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CHAPTER 1 


INTRODUCTION AND SUMMARY 


1.1 INTRODUCTION 

The ever increasing demand on spacecraft power systems for 
improved efficiency and reliability, smaller size and lighter 
weight, coupled with continuous growth in dimension and com- 
plexity of the spacecraft’s payloads, has focused attention 
on a major deficiency - the ability to design, test, and 
trouble-shoot large-scale power systems. 

During the past several decades, numerous efforts have 
been made to develop new and more powerful techniques for 
modeling and analysis of spacecraft power components (equip- 
ments) and subsystems. However, when the components and 
subsystems are interconnected to form a complex system, it 
is quite difficult to predict the total system response even 
though the behavior of individual components may be well un- 
derstood and documented. This is due to many undesired 
interactions that exist among highly nonlinear components. 
The need of a comprehensive power system modeling tool is, 
therefore, most critical since the elaborate design verifi- 
cation through integrated systems hardware testing is pro- 
hibitively expensive or often impossible. This is 
particularly true in the future space station system since, 
due to the large dimension and complexity, the complete sys- 
tem can only be assembled in space. Therefore, a comprehen- 
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sive computer model that can actually predict a system's 
local and global behaviors is most critical for the success 
^of future missions. 


1.2 OBJECTIVES 

The study objectives in Phase 1, 'Large-Signal Modeling 
and Analysis of Spacecraft Power System' are described as 
follows: 

1. Incorporation of all component models generated under the 

previous NAVY contracts into CAMAPPS. 

The large-signal models developed under the previous 

NAVY contracts consist of: 

• Solar array models - an analytical model and an em- 
pirical data model for the solar array dc models. 

• Shunt regulator model - a particular shunt circuit 
model for dc, small-signal and large-signal analysis. 

• Switching converter model - A buck converter with 
constant frequency and single-loop control. 

• Payload models - reactive load and empirical data 
model . 
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2. Generation of large-signal componet models and the macro 
component model library. 

Large-signal component models need to be generated to 
facilitate analysis and simulation of comprehensive 
spacecraft power systems. The component models are in 
the form of macro-models which can be easily configured 
into a particular spacecraft power system. The macro- 
model library to be generated shall include: 

• Solar array 

• Solar array switching unit and control 

• Shunt regulators 

• Battery chargers/dischargers 

• DC-DC converters 

• Payload 

3. Analysis and simulation of a spacecraft power processing 
system 

Large- signal behaviors of solar array power systems 
will be analayzed. A spacecraft power system will be 
configured and simulated using various type of component 
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modules to demonstrate the utility of the proposed 
CAMAPPS . 

4. Transfer of software program of CAMAPPS to NASA 

5. To identify future development effort on CAMAPPS and 
possible integrated system configuration for CAMAPPS 
executive, data base management and user interface. 

1.3 SUMMARY OF THE WORK 

The proposed works for Phase 1 are accomplished and sum- 
marized in the following: 

Task 1. Component model library for switching converters - 
Chapter 3 

Since there exist many different types of power stage 
topologies and control methods, the large-signal behaviors 
of a regulator depend largely on the type of power circuit 
topology and control. Thus, for maximum flexibility, it is 
best to develop models for each functional block as inde- 
pendent modules. A regulator can then be configured by col- 
lecting appropriate pre-defined modules for each functional 
block. Also, any of these block modules can be freely re- 
placed with alternatives, it provides maximum flexibility in 
building a variety of regulators. The macro modules devel- 
oped for each functional block are in the following. 
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Converter topologies 


Buck 


Boost 


Buck/boost (Flyback) 


Forward 


Control methods 

• Single-loop and multi-loop feedback with the following 
duty-cycle controls 

— Constant frequency 


— Constant T 


on 


— Constant T 


off 


• Current-injected and SCM control with 


— Constant frequency 


Constant T control 
off 
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Task 2. Other component model generation 


In order to complete the component model generation for a 
comprehensive spacecraft power system, the following macro- 
modules are developed. The capability of each module is 
demonstrated using a simplified Direct Energy Transfer (DET) 
system. 

• Solar array switching unit and control - Chapter 4 

• Shunt regulators - Chapter 5 

• Battery discharger - Chapter 6 

Task 3 . Analysis and simulation of spacecraft power proc- 
essing system - Chapter 7 

Large- signal behaviors of solar array power systems are 
analyzed. Stability of the solar array system operating 
points with a nonlinear load is analyzed. The state-plane 
analysis illustrates trajectories of the system operating 
point under various operating condition. Stability and 
transient responses of the system operating near the solar 
array's maximum power point are also analyzed. The solar 
array system's mode of operation are described using the Di- 
rect Energy Transfer spacecraft power system. The DET system 
is simulated for various system operating conditions. 
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Task 4. Transfer of software program of CAMAPPS to NASA/GSFC 


Software program 

The software program developed under Phase 1 effort is 
ready to be transferred on tape to NASA 

User's handbook 

The user's handbook is included in this report and con- 
sists of the following documentation 

— Macro-model generation 

— Macro-model description 

— User's input requirement 

— Macro-model output 

— Sample modeling and simulation 

— Program listings 
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CHAPTER 2 


MODELING APPROACH OF LARGE-SCALE SPACECRAFT POWER SYSTEM 
2.1 INTRODUCTION 

Considering the size and complexity of today's spacecraft 
power processing systems, the digital computer is perhaps the 
only viable tool for system modeling and simulation. Two 
general classes of power system modeling programs appear in 
the literatures. These are : 

• Generalized circuit and system analysis models. 

• Dedicated models of specific systems. 

The generalized circuit analysis models, as a class, are 
represented by programs such as SPICE, SCEPTRE, SYSCAP, ICAP, 
ASTAP and a number of others. All are similar in that they 
offer a set of typical device modules which are assembled 
into a network of interconnecting nodes. This class of com- 
puter programs shares a common drawback, such that with an 
increasing number of components, the memory requirements in- 
crease dramatically. By use of dynamic memory management it 
may be possible to analyze a large-scale system, but only at 
the cost of a geometric increase in computation time. The 
labor involved in model-building is rather extensive and the 
resulting model, in general, lacks flexibility for future 
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modifications. Furthermore, the analysis capability of these 
canned programs is limited and usually not effective for 
design aid or trouble-shooting. 

Dedicated models of specific power systems have been used 
successfully but only on a limited basis. Typically this 
class of model is capable only of dc analysis. Execution 
speed and cost of such programs tend to be much less than that 
for generalized circuit and system routines [1]. 

One way to develop a model with such attributes is to 
modularize various components of the power system. There- 
fore, a comprehensive power system model must be capable of 
performing various analyses: dc, small-signal and large- 

signal transient. The model should have the following at- 
tributes : 

• flexibility for future modification of component models 
and system configuration, 

• accuracy to provide sufficient details about the behavior 
of individual components for design and troubleshooting, 

• verifiable whenever equipment or subsystem experimental 
testing is possible, 

• efficiency to minimize computer core memory and computa- 
tion time. 
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2.2 MODELING APPROACH OF LARGE-SCALE SYSTEMS 

As discussed in Chapter 1 , defining a complete spacecraft 
power system entirely in terms of circuit elements is not 
feasible. A practical approach is to consider that a large- 
scale system can be defined by the interconnection of many 
smaller subsystems. The system dynamics, therefore, stem 
from the individual subsystems and their interconnection laws 
that identify the signal flows among them. Due to the com- 
plexity of the individual components and their intercon- 
nection scheme it is usually very difficult to predict the 
system's overall behavior. A desirable analysis or synthesis 
procedure for a large-scale system is the separation of these 
two levels (component and system). One thus prefers an 
analysis or synthesis procedure to be based, on one hand, on 
" local ” subsystems and, on the other hand, on the " global " 
interconnection laws. 

In order to utilize the modular concept, each component 
should, be modeled as an unterminated two-port network because 
the terminal characteristics of both source and load are un- 
known until the complete system is configured. Thus, the 
information flow can be obtained in both directions. The 
interconnection law should not only describe accurately the 
signal flows between the subsystems but also be standardized. 

There exist several kinds of two-port representations such 
as impedance, admittance, chain, and hybrid matrices. The 
impedance or admittance matrices cannot be used as a standard 
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interconnection law because their inputs and outputs are in 
alternate forms. The chain matrix is not as flexible as the 
hybrid matrix because it only characterizes a system in uni- 
direction. 

Figure 2.1 illustrates the interconnection law for unter- 
minated modules. Each component, in general, receives two 
inputs, and the preceding and following modules, 
and provides two outputs, i^ and to the adjacent modules. 
The terminal characteristics of each component are described 
by the matrix equation called "hybrid g-parameters" in 
Eq.(2.1). 


’ ^1 ‘ 


' ^11 

^12 ' 


' ^1 ' 


. ^2 . 


^21 

1 


. . 

( 2 . 1 ) 


For modeling a two-port, one can always define the output and 
input variables as shown in Fig. 2.1. For example, a single 
stage LC low-pass filter can be modeled with state equations 
as shown in Fig. 2. 2. Since these unterminated two-port com- 
ponents have the same interconnection law, they can be 
interconnected arbitrarily to form a system. 

Having identified the two-port component coupling ap- 
proach, efforts have been made to find a software system 
which can be adapted to multi-port component modeling and 
system integration. The program EASY5 developed by Boeing 
Computer Service has been selected as the host software sys- 
tem. 
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Figure 2.1 Two-port subsystem coupling using the hybrid G-matrix 
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Figure 2.2 Example of an unterminated two-port modeling 
with state equations 
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2.3 EASY5 SYSTEM MODEL GENERATION AND ANALYSIS 


For a large-scale spacecraft power system, the large- 
signal performance evaluation is invariably limited to tac- 
tics closely identified with time-domain simulation 
techniques. The actual operating point trajectories or 
transient responses due to the large-signal disturbances can 
be simulated by a digital computer. In this way, the exist- 
ence of unstable solutions can be observed along with peak 
values and actual waveforms. Thus, the model development of 
each nonlinear component of a spacecraft power system in- 
volves two separate models, the large-signal time-domain 
simulation model and the small-signal linearized model. 

For time-domain simulation, the EASY5 numerical inte- 
gration routine calculates current values of the states and 
output variables for each component from the input values. 
It then passes along the output values according to the 
interconnection law. There are seven different types of nu- 
merical integration routines available in EASY5. Each offers 
different numerical stability and computational efficiency 
dependent on the type of nonlinearity of a system. To exam- 
ine the approximate large-signal behavior of a system ana- 
lytically, rather than by time-domain simulation, one may 
simplify the system focusing on a particular component of 
interest. For instance, the large-signal behavior of the 
solar array system's operating point can be analyzed by sim- 
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plifying the system to a second order nonlinear system as 
illustrated in Chapter 7 . 

The EASY5 software [2] provides a modular approach to dy- 
namic system model building and analysis. It was chosen 
particularly for its effective means of assembling complex 
component modules into a system. Predefined component mod- 
ules can be stored in the EASY5 component library. The com- 
ponent modules can have several input and output ports and, 
in each port, information can flow into and out of the port. 
All detailed connections of signal paths between component 
blocks are programmed with minimal user intervention. 

Once all the necessary components are modeled, they are 
stored in the EASY5 macro component library. As shown in 
Fig. 2. 3, the EASY5 consists of two programs. Model Generation 
Program and Analysis Program, and a library of predefined 
Standard Components. Predefined Standard Component models 
include many of the common effects found in dynamic systems 
such as standard analog and discrete controllers (propor- 
tional, integral, differential, etc.) transfer functions, 
nonlinearities (relay, saturation, hysteresis and switches), 
linear state equation model, etc. 

Model generation : The EASY5 Model Generation Program uses 
a block diagram approach for constructing various system 
models. All interconnections between the component models 
are accomplished by the Model Generation Program according 
to the user specified Model Description Data. This is ac- 
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complished by matching the input and output port quantities 
of each interconnection. The program also produces a com- 
plete list of input data that is required by each component 
to complete the model description. The scalar and vector 
parameters and tabular data required for the analysis are 
included in this list. 

Analysis : Once a system model has been constructed, the 
Analysis Program performs various dynamics analyses according 
to the user provided Analysis Data. The Analysis Data in- 
cludes parameter values, initial conditions of states, vari- 
ous analysis controls and commands. Among the various 
analysis controls in EASY5, it is important to note that the 
INTEGRAL control allows one to perform the analysis at vari- 
ous levels such as system, subsystem, or component. The 
analysis in the component or subsystem level is done by 
'freezing' the states of the subsystem of little interest. 
This control feature is particularly useful when one is in- 
terested in investigating behaviors of a component or a sub- 
system in the presence of the entire system. A list of the 
analysis techniques available through the Analysis Program 
are given as follows. 

• Time history generation - nonlinear 
simulation 

• Linear model generation with eigenvalue and 
eigenvector calculations 
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Frequency response analysis 
Root locus calculation 
Steady-state analysis 
Stability margin calculations 
Optimal controller synthesis 



CHAPTER 3 


SWITCHING REGULATOR MODELING 

3.1 INTRODUCTION 

Because of high efficiency, small size and light weight, 
switching regulators are widely used in spacecraft power 
systems. Switching regulators can be characterized, as shown 
in Fig. 3.1, by the three basic functional blocks: power 

stage , analog feedback controller , and digital signal 
processor . 

The power stage consists of energy storage elements and 
switches. Transfer of the input power to the load is con- 
trolled by the duty ratio of the switch (transistor). The 
analog feedback controller usually contains an error ampli- 
fier and a compensation network. For a single-loop control 
it senses the output (load voltage), and for a multiple- loop 
control it senses the output and states (i.e., load voltage, 
inductor current, capacitor voltage). The digital signal 
processor includes a ramp generator, a comparator, and 
latches. It takes the control voltage from the analog feed- 
back and converts it to the pulse-width-modulated 
(PWM) signal which controls the switches in the power stage. 



Figure 3.1 Switching regulator model 
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Since there exist many different types of power stage 
topologies and control methods, the large-signal behaviors 
of a regulator depend largely on the type of power circuit 
topology and control. Thus, for maximum flexibility, it is 
best to develop models for each functional block as inde- 
pendent modules. A regulator can then be configured by col- 
lecting appropriate pre-defined modules for each functional 
block. Also, any of these block modules can be freely re- 
placed with alternatives, providing maximum flexibility in 
building a variety of regulators. 


3.2 POWER STAGE MODELING 


Generally, three different circuit topologies exist in the 
course of a complete switching cycle. A converter is char- 
acterized by three sets of state equations during a switching 
cycle as follows: 

X = A^x + B^u 

y = C^x + E^u , for interval (3.1) 

X = A2X + B2U , 

y = C 2 X + E 2 U , for interval (3.2) 

X = A^x + B^u , 

y = C^x + E^u , for interval T ^2 (3.3) 

Here x is the state vector, u the input vector and T the 

s 

switching period. A^ and B^ (i=l,2,3) are square matrices 
which describe the three circuit topologies and the effects 
of the input vector, u. These three matrix equations can be 
combined into one by using the three switching functions 
shown in Fig. 3. 2. With use of these switching functions, a 
single state equation suffices to describe the converter, as 
in Eq. (3 . 4) . 
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X 


3 

= Z [A.x + B.u]d. (t) , 

i=l ^ 1 ^ ' 

3 

y = I [C.x + E.u]d.(t) , i = 1,2,3 (3.4) 
i=l ^ ^ ^ 

where d^ * s are either zero or one depending on the specific 
time interval involved. 

Equation (3.4) is not only discrete but also nonlinear 
because d^(t) is a function of the state vector, x, and the 
input vector, u. The power stage model shown in Eq.(3.4) is 
directly implemented with the appropriate switching function 
for the time-domain model. The intervals defined above in- 
clude the possible inductor-current discontinuous-conduction 
mode (DCM), as illustrated in Fig. 3. 2(b). For the 

continuous-conduction mode (CCM), there are only two switch- 
ing intervals, T^j^ and Tj,^,as shown in Fig. 3. 2 (a). In the 

time-domain model, these two different modes of operation 
(DCM or CCM) are naturally determined depending on circuit 
parameter values and input and output requirements. 
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3.2.1 EASY5 MODEL FOR CONVERTER POWER STAGES 


As described in the previous section, power stage, in 
general, has three sets of equations. These three sets of 
matrices coefficients are stored in the power stage module 
and appropriate sets of matrices are selected according to 
the status of d(t) which is determined by the digital signal 
processor module. As illustrated in Fig. 3. 3, when d(t) is 
equal to one, (interval the first matrix set, A^, B^, 

C^,and D^, is used to calculate state equations, the second 
matrix set is used when d(t) is zero, (interval and the 

third matrix set when both inductor current and d(t) are 
equal to zero (interval ) • 

The state variables are x = where ij^ is the 

energy- storage inductor current and v^ is the output-filter 
capacitor voltage. For the buck/boost (flyback) converter, 
X = [0 v^] where 0 is the flux in the two-winding energy- 

storage inductor. Input vectors are u = [v^ ^2^' 
is the input voltage and ±2 load current. 

All switches (transistors and diodes) are regarded as 

ideal switches, and the winding resistance (R^^) of the 

inductor and ESR (R ) of the capacitor are included in the 

o 

power stage model. 
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For the buck power stage denoted as [BC], the three, lin- 
ear equivalent-circuit models can be generated according to 


the time intervals and Tp 2 , which are illustrated 
in Figs. 3. 4. The power stage is then described by the 
piecewise linear equations for each subinterval as following: 


For 
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interval T ^2 


Figure 3.4(b) Three switched equivalent circuits 
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For the following converter topologies, EASY5 macro mod- 
ules are also generated. Detailed macro-model descriptions 
are given in the User's Handbook. 


[BT] boost-converter power stage 

[FB] buck/boost (flyback) converter power stage 

[FW] forward-converter power stage 
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3.3 ANALOG FEEDBACK CONTROLLER 


The commonly used analog feedback controller is basically 
linear for normal modes of operation. The controller can be 
modeled by either transfer functions or state equations. 

The transfer- function model provides more design insight 
so both analysis and synthesis of the controller can be done 
quite easily. Pole-zero compensations of the closed-loop 
system can be directly implemented using the Bode plot tech- 
nique. However, the second-order effects of the feedback, 
such as the steady-state switching ripple component, is lost. 

Modeling with state equations can describe the circuit in 
more detail and can easily include the nonlinearities such 
as op-amp saturation and protection circuits. However, the 
circuit topology must be pre-determined and, generally, the 
model requires more memory space. Also, the model does not 
provide direct information about the location of poles and 
zeros to facilitate design of the feedback compensator. 

The choice of model depends on the objective. For design 
and frequency-domain analysis the transfer-function model 
suits better; for more accurate time-domain analysis and 
simulation, the state-equation model is preferable. 

A typical multiple-loop controller for switching regula- 
tors senses the state (inductor current) and output (output 
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voltage). The sum of the outputs from the current- feedback 
loop and the output- voltage loop forms the control voltage 
or error voltage, which is input to the digital signal 
processor. 


3.3.1 EASY5 MODEL FOR ANALOG FEEDBACK CONTROLLER 

The EASY5 model for the analog feedback controller is 
characterized generally by the voltage-feedback module and 
current-feedback module. When a single- loop feedback control 
is used, the current feedback module is omitted. 

3 . 3 . 1 . 1 VOLTAGE -FEEDBACK COMPENSATOR 

The following two commonly used compensators are modeled 
with transfer functions. 

• [MP] two-pole one-zero compensation 

Vj, = (3.5) 

s(l + s/Up) 

v^ : switching regulator output voltage 

K : voltage dividing factor 

E„ : reference voltage of op. amp. 

Vg : compensator output voltage 
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One possible circuit implementation of the above transfer 
function is shown Fig. 3. 5 (a) and its gain characteristic is 
shown in Fig. 3. 5(b). The transfer function of the 
compensator circuit shown in Fig. 3. 5 (a) is 


Vg 1 + sC2(R^+ R 2 ) 

F(s) = — = (3.6) 

Vq sR^C^d + SC2R3) 


By comparing Eqs.(3.5) and (3.6) one can determine the 
circuit parameter values for resistors and capacitors in 
Eq.(3.6) from its design parameter values, “p 

vice versa. For EASY5 modeling, the transfer function in 
Eq.(3.5) is internally converted to state equations. 

) 
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[PZ] 


Two-pole two-zero compensation 


Ve = ^ [Kv^-EeI 

s(l + s/«p) (3.8) 

Figure 3.6 shows one possible circuit implementation of 
the two-pole two-zero compensator and its gain character- 
istic. The transfer function of the circuit is 


(1 + sC^R3)[1 + sC2(R^+ R^)] 

F(s) = ^ (3.9) 

sRiCi(l + SC2R2) 


Internal to EASY5, 
into state equations. 


the transfer functions are converted 





3.3. 1.2 CURRENT- FEEDBACK MODULE 

Two current feedback schemes are modeled; one is referred 
to as standardized control module (SCM), the other is 
current- injected control (CIC). 


• [SM] Standardized control module(S.C.M) ^ ^ 

In SCM, as shown in Fig. 3. 7(a), voltage across the 
energy- storage inductor is sensed through a transformer and i 
is integrated, so that the integrated voltage (Vj) has the 
same shape as the ac component of the inductor current. That 
is, the transfer function from the inductor current to the 
output voltage of the current-feedback loop (Vj) is a simple 
gain. These relations are shown in the following: 

V = n L di,/dt (3.11) 

dC Xj 

Taking the Laplace transform of equation (3.11) 

V (s) = n L s i,(s) (3.12) 

clO u 

Since 

Vi(s) = - v^^(s)/(s C^R^) (3.13) 

from equations (3.12) and (3.13) 

Vj(s)/iL(s) = - n L/(C^R^) (3.14) 
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SM 



Figure 3.7 Current- feedback module 
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As illustrated in Fig. 3. 7(a), the op-amp and capacitor Cl are 
common to both the current loop and the voltage loop. Thus, 
the control voltage (v^) represents the total error voltage 
derived from both the current loop and the voltage loop. 


• [CC] Current-injected control module (C . I . C) 

In CIC, as shown in Fig. 3. 7(b), the switch current is 
sensed through a current transformer. The voltage drop 
across the current-sensing resistor (R^) is proportional to 
the switch current as described in following equations. 


i /n 
s' 



(3.15) 


Thus, the control voltage after summing the output voltages 
of both the current- loop and the voltage- loop becomes 


= - V, 


- V, 


(3.16) 
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3.4 DIGITAL SIGNAL PROCESSOR 


( 


The digital signal processor (DSP ) converts the control 
signal from the analog feedback controller into discrete- 
time pulses to control the ON-OFF of the power switch. The 
DSP includes the threshold detector (comparator), the ramp 
function and the latch circuit. 

The threshold detector compares the control signal, 
and the reference signal (either dc or a ramp) to generate 
the trigger signal for the switch command. For the multiple 
loop control, the control signal, is the sum of the out- 

puts from all feedback loops. In this case, the waveform 
already includes a ramp function (proportional to the 
inductor current or the switch current). The external ramp 
function is optional, it is needed to stabilize the system 
operating at a constant frequency and a duty ratio greater 
than fifty percent. Figure 3.8 illustrates the PWM model for 
a constant-f requency control. For each different duty- 
ratio-control law, such as constant frequency, constant T^j^ 
and constant the PWM model is different. These dif- 

ferent modules are stored in the component module library, 
as in the case of power stage models. 
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3.4.1 EASY5 MODEL FOR DIGITAL SIGNAL PROCESSOR 

The digital signal processor is modeled according to its 
control law such as constant- frequency control-- [WM] , con- 
stant on-time control-- [NN] , and constant off-time 
control-- [ FF] . 

The protection functions are also included in the model, 
such as the peak-current limiter for power components and the 
converter shut-down in the event of sensed abnormalities such 
as over-voltage, under- volt age, or over-current beyond a 
predetermined and tolerable level. Nonlinearities, such as 
op-amp saturation and duty-cycle limiter, are also included 
in this module. 


[WM] constant- frequency control 


For a single-loop control model, as shown Fig. 3. 9(a), an 
external ramp is used to implement the duty-cycle signal and 
the frequency of the external ramp function is identical to 
the switching frequency. For two- loop control, as shown in 
Figs. 3. 9(b) and (c), the inductor-current or switch-current 
information determines the point at which the control voltage 
(Vj^) intersects the reference voltage (Vj^). The reference 
voltage (Vj^) is simply a fixed threshold voltage (v,jjj) when 
an external ramp is not used. If an external ramp is used 
to solve the instability problem, which occurred when the 


41 



duty ratio exceeds fifty percent^ the reference voltage 
becomes the sum of the threshold voltage and ramp voltage. 

The FORTRAN implementation of this module is shown in the 
flow diagram. Fig. 3. 10. It consists of four basic elements: 
the ramp- function generator, selection of the control method, 
comparator and protection. To generate the ramp function, a 
normalized time (Tj^) is defined by using the EASY5 internal 
variable TIME. 

= (TIME + Tg)/ Tg 

where T is the switching period. By subtracting integer N 
which is generated by the library function IDINT(Tjj) from T^j, 
the normalised ramp £unotion(T^ - N) is produced. Then the 
ramp function is generated by multiplying the desired 
amplitude. These waveforms are shown in Fig. 3. 11. 

Since this module can handle the single- loop control and 
multi-loop control (with or without external ramp function), 
one of these options is selected and the appropriate refer- 
ence voltage(Vj^) and control voltage(Vj^) are assigned. The 
reference voltage and the control voltage are compared to 
determine the status of the switch according to the switching 
logic illustrated in Fig. 3. 9. Whenever one switching period 
ends, the switch is reset to the ON position and the process 
is repeated for the next cycle. Duty-cycle limiter and 
overcurrent-protection logic are included in this module. 
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If any aforementioned abnormalities are sensed, the shut-down 
command will turn off the switch. The switch status (ON or 
OFF) determined in these processes, d(t) , is . transferred to 
the power stage to execute the corresponding state equations. 

• [FF] constant-OFF-time control 

In the constant-OFF-time control, regulation is achieved 
by controlling ON- time interval When the two- loop con- 
trol is implemented, as shown in Fig. 3. 12(b), the inter- 
section of the descending control voltage (v^^) with the 
threshold level marks the end of the ON-time interval. For 
the single-loop control, however, additional elements are 
required because the control voltage has only the output 
voltage information. One possible implementation is illus- 
trated in Fig. 3. 12(a). An external ramp function is initi- 
ated at the end of the constant-off-time interval. The 
intersection of the control voltage with the ramp function 
determines the end of the ON-time interval . 

In the EASY5 model, as shown in Fig. 3. 13, the control 
voltage, v^, and the reference voltage, Vj^, are determined 
according to the choice of the control technique, single- loop 
or two-loop control. If the switch is ON, those two volt- 
ages, v^ and Vj^, are compared to determine the instant of 

turn off (T ). Once the switch is turned OFF, it will remain 
c 

for a pre-determined OFF- time interval (Tq^^.) . At the end 



OFF time 


(b) SCM 











of the the switch is turned ON and the on-time instant 

(Tj^) is reset. 

• [NN] constant-ON-time control 

In constant-ON-time control, the output regulation is 
achieved by controlling the OFF- time interval. For two-loop 
control, as shown Fig. 3. 14(b), the intersection of the as- 
cending control voltage with the threshold level marks the 
end of OFF- time interval For single- loop control, as 

shown in Fig. 3. 14(a), the intersection between the control 
voltage and descending ramp signal decides the turn-on in- 
stance. Figure 3.15 shows the flow chart of constant ON-time 
control module. In this module, control logic is basically 
the reverse of the constant-OFF-time control. Comparison 
between v^ and Vj^ occurs during OFF- time interval, and after 
a predetermined ON-time interval, the switch is turned off 
and the OFF- time instant (T^^) is reset. 
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T„v variable 
ON 

. .OFF time 
constant 


ON 


ON 


ON time 


(b) SCM 




NN CONSTANT ON TIME CONTROL 



Figure 3.15 Flow-chart for constant ON-time control 
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3.5 SIMULATION EXAMPLES 

Since a switching regulator is modeled as an unterminated 
component to give the maximum flexibility to configure a 
system, it is necessary to connect a load for the simulation. 
A simple LCR load model [LO] is developed, as shown in 
Fig. 3 . 16 . 

The state equations are 
= <''1 - '' C > / ^ 

''C = <^L - V"> '' 



Figure 3.16 LCR load 
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As described in Section 3.1, a complete converter circuit 
can be configured by arbitrarily choosing a power stage, a 
analog feedback controller and a digital signal processor. 
To demonstrate the versatility of this modeling scheme, the 
following two examples are given: 

• Example 1 : Buck regulator employing SCM feedback and a 
constant-OFF-time duty-cycle control 

Following modules are selected to model a switching regulator 
system as shown in Fig. 3. 17. 

[BC] buck power stage 

/ 

[MP] compensator (two-pole one-zero) 

[SM] SCM (current loop) 

[FF] PWM (constant-OFF-time control) 

( 

[LO] load 

The EASY5 Model Description Data [BUCKCF.MOD] is shown in 
Fig. 3. 18(a). It includes the calling of each module from the 
Macro-Module Library and the interconnections between the 
component modules. This is accomplished by matching the 
port quantities between the input port and the output port. 
According to the user supplied Model Description Data, the 
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\\t — ll. 


Load 
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EASY5 Model-Generation Program generates the regulator cir- 
cuit configuration (Fig. 3 . 18(b) ) and input data requirement 
list (Fig. 3. 19). Figure 3.20 shows the user's input Analysis 
Program that includes the parameter values, initial condi- 
tions of state variables, and various types of analysis com- 
mands. Figures 3.21(a)-(c) show the steady- state regulator 
responses such as the regulator output voltage, the inductor 
current, and the control voltage vs. threshold voltage. 
Transient responses when load resistance is increased ten 
percent at 1ms are shown in Figs.3.22(a)-(b) . 
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««* BUCK LARGE SI GNAL < MODULAR APPROACH) 

-It «• if *4?- «•****#•«■**■«•*«■■»» -K ■«•#** -tHHt #**•«•* * 

** REVISED 2/4/86 

MACRO FILE NAME=MACROS 

♦LIST MACRO COMPONENTS=BC. LO, FF. SM, MP 

MODEL DESCRIPTION 

L0CATI0N=12. DC. INPUTS=LO( IL=I2) . FF < Q=IQ ) 
L0CATI0N=17, LO. INPUTS=BC( V2=V1 ) 

L0CATI0N=34. SM. INPUTS=BC < VL=VL ) 

L0CATI0N=36, MP, INPUTS=BC(V2=V0) 

L0CATI0N=43. FF, INPUTS=SM(VI=VI ), MP< VE=VE). BC< IL=IL) 

END OF MODEL 

PRINT 


Figure 3.18(a) User inpuL model description data 

I BUCKCF . MOD ] 
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♦♦♦♦•ft***** 

* * 

♦ BC ♦ 
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Figure 3.18(b) EASY5 generated model configuration 

[ BUCKCF . MGL ] 
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Figure 3.19 Input data requirement list [BUCKCF.MGL] 



TITLE=BUCK. sen CONST OFF TINE (LARGE SIGNAL) 
♦REVISED 1/24/86 
♦*BC (POWER STAGE) 

INITIAL CONDITIONS 
IL BC=10*VC BC=20 
PARAMETER VALUES 
L BC=250E-6. C BC=400E-6 
RL BC=1E*2. RC BC- lE-2 
VI BC=28 

**VCG LOOP COMPENSATER 
K MP=. 3. ER MP=6 
GA MP=4488. WA MP”628. WB MP-1E5 
** CURRENT LOOP(SCM) 

NV SM=. 08. C SM=7. 275E-9. R SM=3. 44E3 
♦♦ PUN 

TI FF=20E-6. VP FF=0. VTHFF*. 5 
ER FF=6. CICFF=0. SCMFF=1 
T0FFF=10E-6. EPSFF*=lE-5 
ILMFF=15 
** LOAD 

RA L0=2,RB L0=2. 2. TC L0=100E-3 
L L0=lE-6, C LO-IE 6 
INITIAL CONDITIONS 
VE MP=0 

IL L0=10, VC L0=20 
PRINTER PLOTS 
ONLINE PLOTS 
I NT M0DE=4 
DISPLAYl 
V2 BC 
DISPLAY2 
IL BC 

D I SPLAY4 ( OVERPLOT ) 

VR FF, VC FF 
TMAX=10E-3. TINC=5E-7 
PRATE=200, 0UTRATC^20 
*PRATE=100. 0UTRATE=4 
SIMULATE 
X IC — X 

PARAMETER VALUES 

RB L0=2, 2 

TC L0=222E-3 

DISPLAYl 

V2 BC 

DISPLAY2 

IL BC 

D I SPL AY4 ( OVERPLOT ) 

VR FF. VC FF 
TMAX=100E-6. TINC=2E-7 
PRATE=10. 0UTRATE=5 
SIMULATE 
XIC-X 

PARAMETER VALUES 

TC L0=lE-3 

DISPLAYl 

V2 BC 

DISPLAY2 

IL BC 

DISPLAYS 

VL BC. VM SM 

DI SPLA Y4 ( OVERPLOT ) 

VR FF. VC FF 
TMAX=6E-3. TINC= 2E-7 
PRATE=200. 0UTRATL=10 
SIMULATE 


Figure 3.20 User input analysis program 

[ BUCKCF . ANC ) 
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TIME ( xO. 1 ms ) 

Figure 3.21 Buck converter steady- state response 

(SCM, constant OFF-time control) 

(a) 

output voltage V 2 

(b) 

inductor current i^^ 

(c) 

comparator input voltages v~ & v, 







( ms ) 


Figure 3.22 Buck converter step- load response 

(SCM, constant OFF-time control) 

(a) load resistance R 

(b) output voltage V 2 

(c) inductor current i- 






• Example 2 : Buck/boost (Flyback) regulator employing CIC 
and a constant- frequency control 

Following modules are selected to model the circuit shown in 
Fig. 3 . 23 . 

[FB] buck/boost power stage 

[MP] compensator (two-pole one zero) 

[CC] CIC (current-feedback loop) 

[WM] PWM (constant- frequency control) 

[LO] load 

The EASY5 Model description data and the EASY5 generated 
system configuration are shown in Fig. 3. 24(a) and (b), re- 
spectively. The input data requirement list is shown in 
Fig. 3. 25. Figure 3.26 is the user input Analysis Program 

list. Figures 3.27(a)-(c) show the simulation results at 
steady state. Transient responses when load resistance is 
decreased ten percent at 1ms are shown in Figs . 3 . 28( a) - (b) 



** FLYBACK LARGE SIGNAL (MODULAR APPROACH) 

REVISED 12/8/85 
MACRO FILE NAME=MACROS 

♦LIST MACRO COMPONENTS=FB. LO. CC. SM> MPi WM 
MODEL DESCRIPTION 

LOCATION* 12, FB. INPUTS=LO( IL=I2), WM( IQ=QQ) 
L0CATI0N=17, LO. INPUTS=FB ( V2=V1 ) 
L0CATI0N=34. CC, INPUTS*FB (PHI=IL) 
L0CATI0N=3<£), MP. INPUTS=FB ( V2=V0 > 
L0CATI0N=43. UM. INPUTS^^CC < VI*VI > , MP ( VE*VE > 
END OF MODEL 
PRINT 


Figure 3.24(a) [FLYBCK.MOD] 
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Figure 3.24(b) [FLYBCK,MGL] 
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Figure 3.25 [FLYBCK.MGL] 


TITLE=FLYBACK, C1C<LAHGE SIGNAL) 
PREVISED 12/8/85 
*«FB (POWER STAGE) 

INITIAL CONDITIONS 
PHIFB=. 8E-4. VC FD* 5. 2 
PARAMETER VALUES 

RS FB= OlSiRC FB=. 007.LS FB=82. 6E-6 
NS FB=10. C FB=330E-6. RP FB=-. 02 
LP FB=400E-6. NP FB=22. VI FB=23 
**VTG LOOP COMPENSATER 
K MP=1,ER MP=-5. 2 

GA MP«=952. 4. WA ■ tr^'SAB. 8. WB tiP=43270 
** CURRENT LOOP(SCM) 

*NV SM=. 08, C Srt-7. 275E-9, R SM=3. 44E3 
** CURRENT LOOP < CIO 
RW CC=53. 6, NI CC=100 
NP CC=22, LP CC-400E-6 
«« PWM 

TI WM=20E-6. VP W1=0,VQ WM«5. 2 
ER WM=5. 2, CICWM=1, SCMWM=0 
** LOAD 

RA L0=1.69, RB L0=-1.5, TC L0=100E-3 

L L0=1E-6.C LO=^ie-6 

INITIAL CONDITIONS 

VE MP=“2. 54 

IL L0=3. 077. VC L0=5. 2 

PRINTER PLOTS 

ONLINE PLOTS 

I NT M0DE=4 

DISPLAY! 

V2 FB, PHIFB 
DISPLAY2 
IP FB, IS FB 
DISPLAYS 
PHIFB 

D I SPLAY4 ( OVERPLOT > 

VR WM, VT WM 
TMAX=5E-3, TINC=5E-7 
PRATE=500, OUTRATE^IO 
«PRATE=100, 0UTRATE=4 
SIMULATE 

xic-x 

PARAMETER VALUES 
TC L0=222E-3 
DISPLAY! 

V2 FB 

DISPLAY2 

IP FB, IS FB 

D I SPLAY4 < OVERPLOT ) 

VR WM, VT WM 
TMAX=aOE-6, TINC=2E-7 
PRATE=50, 0UTRATE=2 
♦SIMULATE 
X IC— X 

PARAMETER VALUES 
TC L0=!E-3 
DISPLAY! 

V2 FB 
DISPLAY2 

DISPLAY3 

PHIFB 

D I SPL A Y4 ( OVERPLOT ) 

VR WM, VT WM 
TMAX=5E-3. TINC=2E-7 
PRATE=200, 0UTRATE=5 
SIMULATE 


i 

Figure 3.26 IFLYBCK.ANC] 
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TIME( xO.l ms ) 


Figure 3.27 Flyback converter steady-state response 
(CIC, constant frequency control) 

(a) output voltage V 2 

(b) inductor current 

(c) comparator input voltages v^^ & 
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T I ME ( ms ) 

1 Flyback converter step-load response 
(CIC, constant frequency control) 

(a) load resistance R 

(b) output voltage v^ 

(c) flux <t> 






CHAPTER 4 


SOLAR ARRAY MODELING 

4.1 SOLAR ARRAY MODEL 

4.1.1 INTRODUCTION 

Solar arrays, made of silicon solar cells, are the main 
source of electrical power for spacecraft power systems. 
Because the individual silicon cell used on these arrays is 
small and produces very little power, many solar cells are 
connected in series and parallel combination to provide nec- 
essary power to the spacecraft. 

Several approaches have been proposed for modeling a solar 
array [3]. First, it is possible to model an array by an 
interconnection of individual solar cells (individual cell 
model approach). Advantages of this approach are: the pa- 
rameters for the individual cell model are readily obtained 
using established measurement and calculation procedures and 
the effects of parameter variation can be easily included. 
Shadowing effects, as well as cell faults, can easily be in- 
serted into the array. However, the major disadvantage of 
this approach is that an extremely large amount of computa- 
tion time is required which might easily exceed the total 
computer capacity. 

Another approach is to use a single-cell model to simulate 
the entire array (macro model approach). The major advantage 
of this approach is the ability to minimize computer time and 
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memory space needed to analyze array behavior. Parameters 
can be obtained by incorporating the measurement/calculation 
of the individual cells and their interconnection scheme or 
simply by making terminal measurements of the array itself. 
Neither of these methods takes much more effort than would 
be needed for the individual cell model. One disadvantage 
of the macro model approach is a new set of macro model pa- 
rameters needs to be developed for each different array con- 
figuration. This tends to make the model less flexible than 
desired. Furthermore, shadowing and individual cell faults 
cannot be easily included in the macro model without, again, 
a redevelopment of model parameters. 

An approach which uses a combination of the macro model 
and the individual cell model could be very advantageous. 
This approach retains the flexibility of the individual cell 
model approach while reducing the necessary computer time and 
capacity. Any fault that occurs in an individual cell can 
be simulated. Shadowing effects can be accounted for in much 
the same way as would be done in the individual model. Sen- 
sitivity of parameters can also be easily analyzed by varying 
individual cell parameters. The disadvantage of using the 
combination model is the need to determine at least two sets 
of model parameters -- one set for the group of individual 
cells and one set for the group of macro models. 

The choice of modeling approaches depends on the partic- 
ular modeling objective( s) . In the present study, the macro 
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model approach, shown in Fig. 4.1, is used for an overall 
system analysis. With this type of analysis, it is presumed 
that the solar array model need not have great detail. 

SOLAR ARRAY DC MODEL 

For this study, two types of dc models are generated: the 
analytical macro model and the empirical data model. For the 
solar array's analytical macro model, the standard solar cell 
dc model, shown in Fig. 4. 2, is used [4], It provides very 
good first order predictions of the dc behavior of the 
silicon cell. Literature is available for determining its 
parameter values. An analytical expression that describes 
the terminal characteristics of the solar cell model in 
Fig. 4. 2 is 


I 

o 


^g “ ^do t < 2 ( Vq I^.R^ ) / kT } 

^ ^o 'o*^s > / ^sh 


1 ] 

(4.1) 


I = current output of the cell 
o 

I = light-generated current 

I, = reverse saturation current 

do 

Q = magnitude of the electronic charge 

k = Boltzmann constant 

T = absolute temperature 

= output voltage at the cell terminals 

R = series resistance 
s 

R , = shunt resistance 

sh 
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Based on the cell model, a solar array macro model is con- 
structed by assembling solar cells in parallel and series 
connections as shown in Fig. 4.1. In Fig. 4. 3, the resulting 
equivalent circuit model is illustrated where its elements 
are lumped parameters and functions of the number of cells 

in series, n , and the number of parallel strings, n . 

s p 

The analytical expression for the terminal I-V character- 
istic of the solar array model in Fig. 4. 3 is then 

I = n [I - I , { EXP (Q(V / n + I.R /n )/kT)} 

o pg do'^'o's os'p'' ’ 

-(V /n +I.R /n )/R, ] (4.2) 

'o's os'p'-'sh ^ ’ 

From Equation (4.2), the solar array operating conditions are 
set by the load power demand and the shunt regulator' s lim- 
iting voltage. Other more elaborate, analytical dc models 
exist, such as multiple-element models which include second 
order effects [5]. The multiple-element models are used to 
accurately model the effects of distributed series resist- 
ance, complex voltage dependencies and radiation. Due to the 
complexities of the model it is not practical to incorporate 
these second order effects in the system study. 

The second model generated for this study is based on em- 
pirical data about the solar array terminal I-V character- 
istics. The solar array illumination level and temperature 
can be chosen as the independent variables. Based on four- 
dimensional empirical data, both linear and high-order in- 
terpolation techniques are employed to define the array I-V 
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characteristic curve for any arbitrary temperature and illu- 
mination. Figure 4.8, in Section 4.1.4, shows the interpo- 
lated I-V characteristics based on two interpolated empirical 
data curves for two different illumination levels. In these 
data, temperature is assumed constant. Other dependent var- 
iables for the array characteristics can be added provided 
their effects can be accounted for by the measurement. If 
the measurement data is available, this approach is simple 
and accurate . 

4.1.2 MODEL GENERATION [AR] 

Equation (4.2) can be rewritten as following. 

t < EXP ( k ( / n^ + I^.R^ / n^ ) ) } 

o p g do 'o o's os p 

-(V / n + I.R / n )/R,] (4.3) 

'o's os'p^'sh 

where K = Q / ART 

: Nominal temperature ( 301°K ) 

A : Empirical data-fitting constant 

Since the above equation is an implicit function, that is, 

V = f(I ,V ), solutions I and V can be obtained using a 
o ' o o' o o 

numerical iteration method. In this model, the Newton's 
Method was used. The algorithm of Newton's Method is ex- 
pressed as 


^n.l = ^n - ^^^n> / ^ <^n> 


n = 0,1,2 
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To incorporate the solar array model in a complete system, 
the solar array load current, I^, is chosen to be the inde- 
pendent variable (input variable) and the output voltage, V^, 
to be the dependent variable (output variable) as expressed 
in Equation (4.4). 

f(V ) = n [I - I , { EXP (k (V / n + I.R /n^))} 

'o' p g do ' o' o' s o s' p 

- ( V / n + I -R / n ) / r' ] - I (4.4) 

'o's os'p''sh o ' ' 

As shown ih Fig. 4. 4, an LC filter was included in the ac- 
tual solar array circuit model for the EASY5 program model 
[AR] . This filter can represent either the output filter of 
the solar array or the cable impedance. It is also important 
to note that assigning the output of the solar array model 
as a state avoids the implicit loop problem in the EASY5 
program. With the solar array output voltage (V^), obtained 
from the Newton's Method, two state equations can be written 
to find the bus voltage (Vg) of this model. 

L (dlj^ / dt ) = - R-Ig - Vg 

C (dVg / dt ) = 1^ - 

There are several papers (4,51 which describe a method to 
include the temperature effect in the solar array model. 



One simple expression for the temperature effect is as fol- 
lows [ 4] : 

AI = f -AT 
c 

AV = ( f + f • R ) AT 

where AT, AI and AV represent the difference between the 
values at a given temperature and the nominal temperature. 
The temperature coefficients of current and voltage are f^ 
and f^, respectively. The current coefficient, f^, can be 
obtained from the empirical data of the solar array I-V curve 
by comparing the short-circuit currents at a given temper- 
ature and the nominal temperature. The is obtained from 
the open-circuit voltage difference at the two temperatures 

with the known f and R„. 

c s 

The illumination change effect of the solar array I-V 
characteristics is modeled by multiplying an appropriate 
factor to the light-generated current, I , under the assump- 
tion that the illumination change affects only the value of 
I . 

g 

4.1.3 EXAMPLE OF SIMULATION 

Figure 4.5 shows the solar array I-V curves at four tem- 
peratures. For this simulation, the model [AR] shown in 
Fig. 4. 6 is used. When the model was simulated, the 
cable/filter section was ignored to obtain the solar cell 
terminal I-V curve characteristics. 


78 


CURRENT [ X O.IA ] 

0.00 0.36 0.72 1.08 1.44 1.80 



Fig. 4. 5 Simulation Result of Solar Array Model 
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MACRO FILE NAME - MACROS 
DEFINE MACRO - AR 

** SOLAR ARRAY MODEL 


ft 

■fi.==r = = = 

===> VP 
===< I X 



■ftftft 

MACRO 

INPUTS 

= IX 

LL 

R 

L 

C 

1 A 

FC 

FV 

NP 

NS 

ft IX = 

INPUT 

CURRENT 


ft LL = 

ILLUMINAI ION 

LEVEL 


* LLS = SLOPE OF ILLUMINATION CHANGE 

* R.L, C, ; CAULE IMPEDANCE 

* TA = AMBIENT TEMPERATURE 

* FC>FV : TEMPERATURE COEFFICIENTS 

« NP = NO. OF PARALLEL ARRAYS 

* NS = NO. OF SERIAL CELLS 

* R.L, C, : CABLE IMPEDANCE 

MACRO OUTPUTS ^ VO IL VB 

* VO = SOLAR ARRAY OUTPUT VOLTAGE 
•K- VB = BUS VOL T AGE 

MACRO CODE 
MACRO STOP SORT 

SOLAR CELL PARAMETERS ***•«• ' 

«■ 

* RS = INTERNAL SERIES RESISTANCE 

* RSH*= INTERNAL. SHUNT RESISTANCE 

* XIO= REVERSE SATURATION CURRENT 
» TN = NOMINAL TEMPERATURE 

» 0 = ELECTRON CHARGE 

» Xl^ = BOLTZMANN CONSTANT 
» XIG= LIGHT-GENERATED CURRENT 

* VOC=- OPEN CIRCUIT VOLTAGE 
«■ 

RS = . 

RSH =- PDO. 
xiG = . 

XIO = 4. 1869E-1 1 
A =■- . 7*7 
TN = SOI. 

Q 1. 60PE -19 
XK = 1. 381E-23 
VOC = . 5D12 

* XKO = O / ( XK * A * TN > 

XK0=39. 3 

«• 

* ILLUMINATION CHANGE WITH A LINEAR SLOPE OF LLS 

FILL = LL. AR — + l.LSAR — * TIME 

ft EFFECT OF ILLUMINATION CHANGE 


Fig. 4. 6 Solar Array Model [AR.MOD] continued 
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XIG = XIG t FILL - 

rr 

Cl = ( 1. RS / RSH ) 

C2 = NP AR— / < NS AR — . 4i- RSH ) 

C3 ^ -NP AR « XIG 
A1 = XKO / NS AR — 

A2 = XKO * RS / NP AR — 

IF ( TIME . NE. 0. ) GOTO +i-HO 

t> 

r. INITIAL GUL-.SS FOR SOLAR ARRAY OUTPUT V0L1 AGE 
VO AR— vn AR — 

li- *** -R- *•»•*•»■*«■•»•»•**■«■«■«•«■ -a- it it -a- 

R-it-tt NEWTON ITERATION 

R * R « R -tt it R R R R R R it -tt R it it * it R it it it R R «• it R R it •«■ it it R * it it * R R R it «■ R R R R it R R 

t-Ri-10 CONTINUE 

VDP = VO AR — 

FV = Cl i;- IL AR — + C2 R VO AR — + C3 > NP AR — 

8< R XIO R DEXP < Ai R VO AR — + A2 R IL AR — ) 

DFV = C? + Ai R NP AR — R XIO R DEXP ( Al R VO AR — 

?< t A2 R IL AR — ) 

VO AR — = VO AR — - FV / DFV 
ZZ = < VO AR — - VOP ) / VO AR — 

IF < DARG ( ZZ ) . LE. l.E-'l ) GOTO ++ 1-?0 

GO TO +t t-10 

RRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRR 

R TEMPERATURE CORRECTION FOR S. A. I-V CURVE R 

R DEL(I> = FC R DEL(T> R 

R DEL(V) =■• < PV + FC R RS ) R DEL(T) R 

RRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRR 

+++20 DELI - FC AR — R < TA AR — - TN ) 

DELV = ( FV AR — + FC AR — R RS ) R ( TA AR TN ) 

ILT = IL AR — + DFLI 
VO AR — -- VO AR — + DFL V 

R 

VL = VO AR — - IL1 R R AR — - VB AR — 

MACRO DERIVATIVES. IL AR — = VL / L AR — 

MACRO DERIVATIVES, VB AR — = ( IL AR — - IX AR — ) / C AR - 
MACRO RESUME SOR'I 
END OF MACRO 

itititititRititRititRRRRitRRititRititRititRititititititRitRititRRRRRRRRRRRRRRRRRRRR 

MODEL DESCRIPTION 
LOCATION = 22. AR 
END OF MODEL 
PRINT 


Fig. 4. 6 Solar Array Model [AR.MOD] 
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This simulation result was similar to the original graph 
drawn from the experimental data in [6]. The solar cell 
model's parameter values were obtained from the experimental 
datja of a 2 X 2cm, lOohm-cm, 8mil silicon cell [4]. The 
analysis program listing [AR.ANC], including all parameter 
values, is given in Fig. 4. 7. 

4.1.4 MODEL GENERATION WITH EMPIRICAL DATA [AE] 

The empirical data model can be generated using the EASY5 
function generator facility. The EASY5 function generator 
constructs a continuous function from a given data set using 
the interpolation technique. The function generator can in- 
terpolate a set of data which is a function of maximum of 6 
independent variables. For instance, the solar array output 
voltage is chosen to be the dependent variable and the solar 
array output current to be the independent variable. A con- 
tinuous I-V curve function can be obtained based on the input 
data set of V vs. I. 

To include the illumination and temperature effect in the 
solar array model, three-independent variable function gen- 
erator, that is, EASY5 standard component [F3], is needed. 

Figure 4.8 shows an interpolated I-V curve created from 
two different data sets belonging to two illumination levels 
of a solar array. In this case, since there are two inde- 
pendent variables, such as the solar array voltage and illu- 
mination level, the function generator, [FV], was used. 
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1 1 TLE^NAS A ( LARGE 3 I GNAL ) 

»T1TLE SOLAR ARRAY 
PARAMETER VALUES 

R AR .001, C AR = 5E-4, L AR =- lOE-6 
NP AR = 1, NS AR=1, TA AR == P.5i. LL AR 1 

PC AR = 8. E-5, FV AR = -2. E-3 

■J*^^#** 

INITIAL CONDI! IONS 
VB AR = O. D 

PRINTER PLOTS 
ONLINE PLOTS 
I NT MODE =--4 
SI MANUAL SCALES 
D1 SPLAY 1 


ILIAR, 

VS, 

VOIAR, 

XRANGE = 

0, 

8, 

YRANGE 0, 0. 16 

IL AR, 

VS, 

VO AR, 

XRANGE = 

0, 

, 8, 

YRANGE = 0, 0. 16 

IL2AR, 

vs. 

V02AR, 

XRANGE - 

0, 

8, 

YRANGE = 0, 0. 16 

IL3AR, 

VS, 

V03AR, 

XRANGE = 

0, 

. 8, 

YRANGE 0,0. 3 6 

TMAX--^ 

0. 16 

, TINC 

0. 0005 




PRATE 

=-100, 

OUT RATE 

>1 





H0UTRATE2=1, PRATE2=--J 
«Pr<INT2 FROM, , OOD, TO. . 006 
SIMULATE 


Fig. 4. 7 User's Parameter Values and Analysis Commands [AR.ANC] 

* 
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VO(S.P) 


Fig. 4. 8 Simulation Result of Solar Array's Empirical Data Model 


84 




In this model, EASY5 command 'MACRO TABLES' is used to 
form a tabular function. The interpolated output from the 
tabular function can be obtained using the following FORTRAN 
command [2], 'CALL FV (TABNAM, S2, SI, S3, ANX, ANY) ' , where 

TABNAM : name of table 

S2 : output (e.g. solar array current) 

51 : input (e.g. solar array voltage) 

52 : input (e.g. illumination level) 

ANX: degree of interpolation for first independent variable 
(e.g. 1.0) 

ANY: degree of interpolation for second independent variable 
(e.g. 1.0) 

The program listings, [AE.MOD] and [AE.ANC], for this 
simulation (Fig. 4. 8) are shown in Figs. 4. 9-10. 
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MACRO FILE NAME = MACROS 

* SOLAR ARRAY MODEL FROM EMPIRICAL DATA 

DEFINE MACRO = AE 

MACRO TABLES TAE(250)2 
MACRO INPUTS ILL 

* ILL = ILLUMINATION LEVEL 

MACRO OUTPUTS 131 IS2 I S3 

DUM VS 

* ISl = SOLAR ARRAY CURRENT AT ILLUMINATION LEVEL 

* OF 1 . O 

* IS2 = SOLAR ARRAY CURRENT AT ILLUMINATION LEVEL 

* OF O. U 

* IS3 = SOLAR ARRAY CURRENT AT ILLUMINATION LEVEL 

* OF INPUT VALUE. ILL. 

* DUM = DUMMY STATE FOR AVOIDING ERROR MESSAGE 

* FROM MODEL GENERATION PROGRAM 

MACRO CODE 
MACRO STOP SORT 
XIN = TIME 


CALL FV 
ISIAE — 

< 

TAEAE — , 
OUT 

OUT, 

XIN, 

1. . 1. . 

1. ) 

CALL FV 
IS2AE — 

7 

TAEAE — , 

our 

OUT. 

XIN. 

0. 8. 1. . 

1. > 

CALL FV 
IS3AE — 
VS AE — 

7 

TAEAE — . 

OUT 

TIME 

OUT, 

XIN. 

ILLAE — . 

1. , 


MACRO DERIVATIVES, DUMAE — = VS AE — 

MACRO RESUME SORT 

END OF MACRO 

MODEL DESCRIPTION 

LOCATION = 33, AE 

END OF MODEL 

PRINT 


Fig. 4. 9 Solar Array Empirical Data Model [AE.MOD] 
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T1TLE=EMPIRICAL DATA MODEL FOR SOLAR ARRAY 
PARAMETER VALUES 
ILLAE = O. 9 

TABLE, TAEAE. AA. 2 

* SECOND INDEPENDENT VARIABLE ( ILLUMINATION ) 

O. 8, 1. O 

* FIRST INDEPENDENT VARIABLE ( VOLTAGE ) 

0. 0. 1. O. 5. 6, V. O, 10. . 

11.5, 14.5, 15.67 ,18.01, 19.85, 

20.65, 21.01, 22.1, 23.01, 24.06, 

25. 02. 25. 52, 26. 02, 26. 34, 26. 7, 

27.0, 27.2, 27.41, 27.6, 27.81, 

28. 0, 28. 2, 28. 4, 28. 6, 28. 8, 

29.02, 29.21, 29.4, 29.61, 29.81, 

30.01, 30.33, 30.54, 30.71, 31.01, 

31.15, 31.33, 31.50, 31.63 

* DEPENDENT VARIABLE < CURRENT ) 


36. 55 

, 36. 5, 

36. 4 

, 36. 35 

i, 36. 3, 

36. 29 

, 36. 2, 

36. 15, 36. 1 

, 36. 02 

35. 95 

, 35. 9, 

35. 7 

, 35. 4, 

34. 8, 

33. a. 

33. 0, 

32. 0, 

31. 2, 

30. 2. 

29. 2, 

28. 5, 

27. 7, 

26. 9, 

26. , 

25. 1, 

24. 1, 

23. 1, 

22. 0, 

20. 9, 

19. 6, 

18. 4, 

17. 2, 

15. 8, 

14. 4, 

13. 0, 

10. 6, 

9. 0, 

7. 7, 5. 

3, 

4. 1, 

2. 6, 1. 

2, O. 

O, 


45. 5, 

45. 5, 

45. 5, 

45. 5, 

45. 41, 

45. 4, 

45. 31, 

45. 3 

, 45. 21 

, 45. 1, 

45. 0, 

44. 9, 

44. 6, 

44. 2, 

43. 3, 

41. 8, 

40. 7, 

39. 3, 

38. 2, 

36. 9, 

36. 5, 

34. 8, 

33. 8, 

32. 8, 

31. 7, 

30. 6, 

29. 5, 

20. 3, 

27. 0, 

25. 7, 

24. 3, 

22. 9, 

21. 5, 

20. 0, 

18. 5, 

16. 9, 

14. 3, 

12. 6, 

11.2, 

8. 6, 

7. 4, 

5. 8, 4. 

3, 3. 

1 



PRINTER PLOTS 
ONLINE PLOTS 
I NT MODE = 3 
DISPLAYl <OVERPLOT) 

ISIAE, VS, VS AE, XRANGE=0, 31. 63 
IS2AE, VS, VS AE, XRANGE=0, 31. 9 
IS3AE, VS, VS AE, XRANGE=0. 31. 9 
TMAX=31. 6, TINC=0. 2, PRATE=10 
SI MANUAL SCALES 
OUTRATE = 1 
SIMULATE 


Fig. 4. 10 User's Parameter Values and Analysis Commands 


[ AE . ANC ] 
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4.2 SOLAR ARRAY SWITCHING UNIT MODEL 


4.2.1 INTRODUCTION 

The solar array switching regulation has been proposed as 
a replacement for conventional regulation techniques to re- 
duce weight, thermal dissipation, and power system complexity 
for large power spacecraft and advanced space missions. 

The solar array switching unit (SASU) contains a voltage 
and current regulating feedback loop. The SASU differs from 
the switching regulator in that it controls the bus voltage 
and battery charge current by switching individual solar ar- 
ray segments as required by the load profile and the state 
of charge of the battery. 

There are two major design parameters for SASU control. 
One is the size of the array segment switching at each in- 
stant. The size of the array segments must be properly cho- 
sen to meet dc condition requirements (e.g. voltage 
regulation specification) of the system. The other design 
parameter is the frequency of switching. The SASU's ability 
to respond to perturbations is a function of the switching 
frequency and the size of the array segment. 

There are four basic configurations (Fig. 4. 11) for the 
SASU and advantages and disadvantages of each are discussed 
in [7]. In this study, the series-switching, parallel-array 
configuration was selected for its low stress on the solar 
array during shadowing, low wiring complexity and simple 
switch drive electronics [7]. 
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Based on the series-switching, parallel-array configura- 
tion, the SASU and control scheme can be constructed as shown 
in Fig. 4. 12. The solar array strings, SSI through SS14, are 
switched-array sections. There is also an active array sec- 
tion which is connected to the main bus without a switch. 
This active array section performs fine regulation of the bus 
voltage while the switched-array sections give coarse regu- 
lation. 

The combined operation of the active and switched-array 
sections is described by the operation of the system during 
exit from an eclipse. When the spacecraft is in an eclipse, 
the array provides no output and the bus voltage is less than 
the setting of the shunt regulator. Given these conditions, 
the shunt is turned off and the current sensors are producing 
a 'count up' signal so that all of the array are switched on. 

As the spacecraft exits eclipse, the full array is on and 
eventually the output current is sufficient to charge the 
battery and supply the load demand. The bus voltage rises 
until it reaches the voltage reference of the shunt regulator 
and the shunt regulator activates to absorb excess array ca- 
pacity, When the shunt current increases to a preset value 
(e.g. 2A), the 'count up' signal is turned off. As the cur- 
rent increases to another preset value (e.g. 5A) the shunt 
current sensors initiate a 'count down* signal which switches 
off the array sections, one by one, until the shunt current 
is stabilized between 2A and 5A. At this time, the switching 
action stops and the shunt regulator regulates the bus with 
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Fig. 4. 12 Solar Array Switching and Control Configuration 






the active array section. Figure 4.13 shows the array char- 
acteristic under static conditions. 

Circuit operation and model performance are illustrated 
using the simulation results shown in Fig. 4. 16 in Section 
4.2.3. 

4.2.2 MODEL GENERATION [AS] 


To generate the solar array switching model, the previous 
solar array model [AR] was modified to add the array switch- 
ing effect. Since the number of parallel solar array strings 
may vary in this model [AS], the solar array I-V curve varies 
accordingly. The array switching depends on the shunt cur- 
rent level and delay time between switchings. Detailed 
switching logic of the SASU system is shown in Fig. 4. 15. The 
additional capacitor. Cl (Fig. 4. 14) , is included to absorb an 
abrupt step change of the solar array switching dynamics. 
The parameter values of R1 and Cl must be properly selected 
to maintain system stability. 

The following state equations are used to implement the 
cable/filter section of the SASU circuit model shown in 
Fig. 4. 14. 


Vl = Vo - R-Il - Vb 

Cl ( dV^^ / dt ) = 

L ( dlj^ / dt ) = 

C ( dVg / dt ) = II - Ix 


where I = I , + I, , 

X sh load 
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: Shunt Current 


NP ; No. of Parallel S.A. 

DL : Delay Time Between Switchings 


Fig. 4.15 Solar Array Switching Logic 
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4.2.3 EXAMPLE OF SIMULATION 


In this simulation, a simple constant-power load model 
[PT] was used for load variation. A full shunt model [FS] 
was included to sense the shunt current. 

Initially 324 parallel solar array strings, each consist- 
ing of 58 solar cells in series, are connected to the main 
bus. They produce a maximum power of 853 watts. A constant 
power (500W, in this example) load is applied to the bus. 
The excess power of the solar arrays will be consumed in the 
shunt regulator. As the shunt current becomes greater than 
5A, the solar arrays are switched off, 20 arrays at a time, 
with time delay of 1ms. The sharp transients of the bus 
voltage are observed at the switch-off and the bus voltage 
approaches the preset value of 28.14V (Fig. 4. 16). For each 
switching step, the shunt current drops by 1.88A, which cor- 
responds to the current of 20 arrays, until it becomes less 
than 5A. Since the bus voltage is regulated, the load 
current remains almost constant. The solar array current, 
Ig^, is the sum of the load current, Ij^, and the shunt cur- 
rent, 

When a system with this model [AS] is simulated, numerical 
integration methods involving variable step- size algorithms 
should be avoided. Since this model includes a routine for 
sensing the present time and comparing it with the delay 
time, the above integration method can not work properly. 
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Fig. 4. 16 Simulation Results of Solar Array Switching System 
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The fixed-step integration methods, such as mode 3 or 4 in 
the EASY5 User's Guide [2], are recommended for simulation 
of this model. 

The system model description program [PSS2 .MOD] , schematic 
diagram, input data requirements list and analysis program 
[PSS2.ANC] are shown in Figs . 4. 17-20 . 


98 


MACRO FILE NAME =• MACROS 


* SASU SYSTEM WITH SHUNT REGULATOR AND * 

» CONSTANT POWER LOAD * 

« REVISED ON 3/3/86 


LIST MACRO COMPONENTS = 
MODEL DESCRIPTION 
LOCATION=l, AS. INPUTS 
L0CATI0N=23. FS. INPUTS 
LOCATION=5. PT. INPUTS 
L0CATI0N=41 . MC. INPUTS 
END OF MODEL 
PRINT 


AS, FS, PT 

= FS<IH=IH). 
= AS(VB=VB) 

= AS(VB=VL> 

= FS<IH=S. 1 ) 


MC(S. 2=IX) 
PT< IL=S. 4) 


Fig. 4. 17 System Model Generation Using SASU [PSS2.MOD] 
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Fig. 4. 18 Schematic Diagram of [PSS2.MOD] 
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PARAMETERS REQUIRED 


cn 

cn 

Hh* 

< 

u. 

Q.Q. 


I 

in 

w 

K 

CC3 

O 

> 

>CL 


cn 

q: 

tu 

H 

UJ 

z: 

< 

q: 

< 

CL 


CL 

Z 


u 

LU 

> 

q: 

o 


< 

o 

z 

o 

cn 

z 

lU 

z 


Q 

z 

< 


(n 

cn 

l-V- 

U 


<£ 

IL 

CUQ. 

£ 




ru 




£ 

3D 

•T 


X 

O 

Q.H 

O 



Q 

Ul 

q: 

D 

O 


U) 

UJ 

H 

< 

H- 

tn 



>< 





UJ 








cr 

X 


o: 








h- 





cn 

< 


< 





A 

£ 


£ 





O 




cn 

cn 


u 

CE 

Q 


a 

< 

u. 

0.0. 

£ 

1“ 

Z 

cn 

z 





Z 

<1 


< 

cn 

n 

33 

n 

o 




Z 

X 

OO 

u 

o 

cr 


a: 






o 

cn 

Q 





a: 

H 


o o 

c: UJ 
<r > 
ui 

a a 
z u. 
< 

< 

cn K 
z < 

Q Q 


cncn 

cn 

HK 

<J 



cn 

<< 

u. 

OO 

c 

K 

z 

< 

o 





o 


z- 

ru 

Ufa 

ra 

Q 


tt 

zo 

cr 


o 

Z 

cn 

> 


O Z 
U UJ 

r 

-j 

< Q 

H Q 
Z 

Z < 


UJ 

£ 

< 

Z 


q: 

UJ 


UJ 

£ 

< 

q: 

< 

0. 


< 


cncn 

cncn 

HH- 

o 


z 

cn 

<< 

OIL 

OO 

£ 

to 


< 


£ 

n 


UJ 

Ui 


A 


33 


H 

H 

A 

GO 

q:»-» 

(no 

O 

< 

< 



cn cn 


cncn 

cncn 


o 

cn 

« 

U.O. 

<c 

OO 

£ 

< 

o 

crcn 


n 

o 

Z-J 

>> 

OH 

cn 

> 


H 





H 

Z 





Z 

UJ 





Ul 

z 





z 

o 





o 

o 

cn 

cn 

H 

(J 

o 

£ 

< 

u. 

O 

£ 

£ 

o 





o 

o 





o 


100 


Fig. 4. 19 Input Data Requirements List 


TITUE-~NASA<LARGt SieNflL) 


SOLAR ARRAY SWITCMING UNIT SYSTEM WITH SHUNT * 

* REGULATOR AND CONSTANT POWER LOAD » 

PARAMETER VALUES 

R AS = .001, C AS = 5E-5, L AS = lE-6, Cl AS = 1 . E-3 
R 1 AS = 1.0 

NP AS = 32A. NS AS=58, NNPAS =20. DL AS = . OOl 

■B-ii- * -fr * <i- *■«•*#* «•«■<{■*«•***«■ « 

INITIAL GONDII IONS 
VB AS = 28. 14 

IL AS = 29. D6S 

VCIAS = 28. 140 
■«•**■«•**■»•«• 

PARAMETER VAL UES 

VR FS --- 28. 14, VSAFS = 15. 

R1 FS = 12400, R2 FS = 2200000, R3 FS = 1.78419 

OM FS = 4. . VTHFS = 3. , ILMFS = 4 5. 

•«■•»•**«•■«•** 

PARAMETER VALUES 
VR PT=15 
SW PT= 20000 
PW1PT=694, PW2P1-500 
TC PT=5E-1,PC PT=2 

Cl MC=1,C2 MC=0. C3 MC=1,C4 MC=0 

B 

PRINTER PLOTS 
ONI.iNE PLOTS 
I NT M0DE=4 
DISPLAY I (OVERPLDT) 

VB AS 

DISPLAY2 (OVERPLOT) 

IH FS, IL PT, 10 AS 
DISPLAYS, NPPAS 
TMAX = 20E-3, TlNC=lE-6 
PRATE = 200, OUTRATE - 30 . 

SIMULATE 


Fig. 4. 20 User's Parameter Values and Analysis Commands [PSS2.ANC] 
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CHAPTER 5 


SHUNT REGULATOR MODELING 

5.1 FULL SHUNT MODEL 

5.1.1 INTRODUCTION 

Shunt regulators are used to clamp the power bus voltage 
when excess solar array power is available. There exist bas- 
ically two types of shunt regulator circuits, linear type and 
PWM type. Modeling techniques for the PWM type have received 
a great deal of attention in the last two decades. A dis- 
cussion of the modeling and analysis of the PWM regulator is 
found in the chapter on The Switching Regulator Model. The 
shunt regulator configurations can also be categorized by 
their functions: full shunt and partial shunt. In a full 

shunt configuration, the shunt elements are placed in paral- 
lel with the load. In a partial shunt configuration, the 
shunt elements only shunt part of the solar array cells. 
Each solar array segment is tapped by a shunt element. In 
this study, the full shunt and partial shunt configurations 
are modeled in detail. 

Figure 5.1 illustrates the simplified block diagram of a 
typical shunt regulator. The basic circuit functions are as 
follows. A divided-down sensed bus voltage is compared to a 
reference voltage. The error voltage is then amplified to 
control the current through the shunt elements. A suitable 
resistor (dependent on the shunt's maximum design current) 
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is used in the shunt current path to monitor current flowing 
through the shunt elements. 

To illustrate shunt regulator modeling techniques, typical 
shunt regulators are modeled in detail in the following 
sections. The circuit diagram and the block diagram are 
shown in Figs. 5. 2 and 5.3. These models take into account 
salient dynamics and nonlinearities including Op-Amp feed- 
back, Op-Amp saturation and darlington cut-off. All diodes, 
including the transistor base-to-emitter junctions, have an 
infinite reverse impedance. The first order model is used 
for transistors. 

Most shunt regulators are composed of two major function 
blocks, the error amplifier (EA) block and the shunt element 
block. 

In some applications, the EA block is implemented with 
three identical amplifier circuits and a majority voting 
.circuit for redundancy. Also, the shunt element block is 
mostly composed of parallel, identical, shunt transistor 
circuits or sequentially turn-on shunt transistor circuits. 
Since it is not practical to generate a generalized model 
which fits any shunt regulator circuit, two typical shunt 
circuits (Type-1 and Type-2) are modeled to demonstrate mod- 
eling procedures. These two circuits represent the most 
commonly used linear shunt regulators in spacecraft power 
systems . 
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BUS 
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Fig. 5.2 Shunt Regulator Circuit Diagram (Type 1) 


bus 
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5.1.2 MODEL GENERATION FOR A TYPE-1 SHUNT 


A. MODELING WITH TRANSFER FUNCTIONS 

A shunt regulator circuit. Fig. 5. 2, is modeled through the 
block diagram shown in Fig. 5. 3. The block diagram shows the 
transfer function representation of each functional block of 
the circuit. This transfer function model is useful when 
creating a generalized model to handle different circuits 
belonging to the same type shunt. The transfer function 
model is particularly useful for the frequency response 
analysis. By providing the values of gains, poles, and ze- 
ros, the user can implement the corresponding models for 
different circuits. 

Although the model is implemented based on transfer func- 
tions, the transfer functions are programmed with state 
equations in the model when the functional block diagram is 
actually modeled. The A^ block, in Fig. 5. 3, can be expressed 
in state equations as 


S2 = Al-Sl + XI 

X°1 = Al-Sl ZO - S2-P0 

The above state variable equations were converted from a 

transfer function which had one pole and one zero. The input 

(SI) of this block is an error voltage between the Kl, V(bus) 

product and V^, and the output voltage is the sum of 

the dc voltage component and the ac voltage component (S2) 

of the A block, 
s 
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SI 


V 


ol 


KI.V 3 - 

S2 + V 

r 


The dununy state variable, XI, which is not a capacitor volt- 
age, is introduced by converting the transfer function into 
state equations. Similarly, the G1 and G2 blocks are ex- 
pressed in state equations. However, in this case, the in- 
puts (SI, S3) of these blocks are and V(bus), 
respectively, and the outputs (S2,S4) are directly fed to the 
summer in order to obtain V ^2 • '^^e LI block represents the 
Op-Amp saturation. This limiting function can be easily im- 
plemented by using FORTRAN 'IF' statements (Refer to program 
listing). The V^, is needed to include the base-emitter 
junction voltage drop of the darlington circuit. The L2 
block represents the limiting of control voltage, V^, to a 
positive value. The K2 block is a gain of control voltage 
to base current of the transistor driver circuit and K3 is a 
factor from the number of multiple shunt elements. The L3 
block represents the vs. characteristic of the shunt 
transistor. The F block represents the functional expression 
of the shunt current, with respect to base current, Ij^. 


'sh = " ^2 

(5.1) 

(5.2) 

where 31 and 32 are the current gains of transistors Q1 and 
Q2, respectively, and of Q1 and Q2 is assumed as 0.7V. 

B. MODELING WITH STATE EQUATIONS 

The same shunt regulator circuit. Fig. 5. 2, is modeled di- 
rectly from the state equation representation of the circuit 
without transfer function blocks. This is a more convenient 
way of modeling a specific circuit since it is not necessary 
to calculate transfer functions. The state variable repres- 
ents the actual capacitor voltage or inductor current. This 
circuit can be divided into three function blocks, i.e., er- 
ror amplifier circuit, error amplifier output circuit and 
shunt element circuit, as shown in Fig.5.4a-c. 

1) Error amplifier circuit 

If the Op-Amp is assumed to be ideal, then the output voltage 
(V^j^) can be expressed in terms of the inputs (Kl*Vg and 
and the state variable, by the following equation. 


31-32-R1-R2 

31-32«R2-R3 + 32-R1-R3 + R1-R2 

0.7 32 (R1 + 31 R2) 

31-32-R2-R3 + 32-Rl'R3 + R1-R2 


- ^ ( K 

R1 ^ ^ 


• V 

1 


) ♦ ( 1 * If )-v, 


The state variable, is expressed as 
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'^01 - ''r - '^cl 
C1*R2 


since the current through capacitor Cl is 


^ol - ^r - '^cl 


In order to take into account the saturation effect of Op- 
Amp, the output voltage (V^^) is limited to a voltage between 
the supply voltage and the ground. 

2) Error amplifier output circuit 

Since the error amplifier circuit, shown in Fig. 5. 4b, is a 
simple linear circuit, the following output voltage equation 
and capacitor-voltage state equation are easily obtained to 
model this circuit. 


V = V - V 
o2 ''b ^c2 


R4 R3 


( R3 + R4 ).V 
J ' c 

R3 R4 


3) Shunt element circuit 


This circuit is basically a shunt transistor circuit with a 

base-drive circuit. If the upper darlington is not in cut- 

» 

off, then the emitter current, I, , to transistor Q1 is 


^b = 


''b - ^ "^02 ^ i 
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where 1.4V is a darlington base-emitter voltage drop and N 
represents the number of multiple shunt elements. The shunt 
current flows only when Vg, which is equal to Vg - (V^2 

1.4V), is positive. Also, before there is any shunt current, 
the voltage across the IKQ resistor must be at least 1.4V or 
must be greater than 1.4mA. Once is determined, the 
shunt current can be expressed in terms of through manip- 
ulation of transistor circuit equations. As a result, 

where Cl and C2 are described in Equations (5.1) and (5.2). 
Finally, the shunt current is limited to the maximum current 
capacity of the solar array current. 

5.1.3 MODEL GENERATION FOR A TYPE-2 SHUNT 

Another typical shunt regulator circuit uses a FET instead 
of a BJT as the shunt element. The following procedures use 
the state equation method for modeling the circuit shown in 
Fig. 5.5. 
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From Fig. 5.5, two equations can be written assuming the 
operational amplifier is ideal. 


- ^r _ ^1 - ^r 


R1 


R2 


V - V 

Vi = Vg - ( gm.Vg^ -H ).R3 


(5.3) 

(5.4) 


From Equation (5.4), can be expressed in terms of Vg 
and resulting in an explicit expression for with an 
assumption of R2 >> R3. 

Vi = Vg - gm.Vg^.R3 (5.5) 
By substituting Equation (5.5) into Equation (5.3), 


gm«V 


gs 


_ ( R1 + R2 ) 
R1*R3 


(V - V 
^ B r 


(5.6) 


The power MOSFET is modeled on the assumption that this de- 
vice is operated in a linear region as shown in Fig. 5. 6. The 
drain current is expressed as 


'd = 


where I, 


gm.(V - V ^) 
^ ' gs gst' 

= drain current 


gm = forward transconductance of MOSFET 


= gate threshold voltage 


If this MOSFET characteristic is included in Equation (5.5), 
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V 

gs 


( R1 R2 ) . 
gm-Rl*R3 ' B 


V 


r 


+ 


V 


gst 


The saturation effect of the Op-Amp is implemented by limit- 
ing the maximum output voltage of the Op-Amp using a FORTRAN 
'IF' statement. The shunt current should be limited from 
zero to the maximum solar array's output current. 


5.1.4 EXAMPLE OF SIMULATION 

For simulation of the Type-1 shunt regulator, system 
[SIMLl.MOD] (FIG.5.8) is configured as shown in Fig. 5. 9. The 
system consists of the transfer-function shunt model [HA], a 
solar array model [AR] and a constant power model [PT]. 

In this simulation (Fig. 5. 7), a constant-power load of 
800W is applied to the solar array which has an output power 
capacity of 853W. After the system reaches a steady state, 
the load step changes from 800W to 750W at 2ms. The decrease 
in the load power causes the shunt current to rise. The shunt 
current is added to the load current so the solar array out- 
put current remains almost constant after a small transient. 
As a result, the bus voltage recovers its specified voltage 
(28.14V) after a transient period. The input data require- 
ments list and analysis program are given in Figs . 5 . 10-11 . 

The same system, using the state-equation model [HB] for 
the Type-1 shunt regulator, is also simulated. 
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Fig. 5. 7 Simulation Results of Type-1 Shunt Regulator Model 
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MACRO FILE NAME = MACROS 

*«-«-***-*«-«*«"«-fr«-»*«-**«-**«-****«-»{- ■»•*«•*«■*■«•■«■«•*«■ «■*■****** 

* SOLAR ARRAY SYSTEM WITH TYPE 1 FULL SHUNT « 

* REGULATOR AND CONSTANT POWER LOAD * 


LIST MACRO COMPONENTS 
MODEL DESCRIPTION 
LOCATION=l, ARi INPUTS 
L0CATI0N=23, HA. INPUTS 
LOCATION=5. PT. INPUTS 
LOCATION=‘<n, MCi INPUTS 
END OF MODEL 
PR 3 NT 


= AR. HA, PT 

= MC(S, 2=IX) 

= AR(VB=VB> 

= AR(VB=VL) 

= HA( IH=^S. 1 ), PT( IL=S, 4) 


Fig. 5.8 System Model Generation Using Type-1 Shunt [SIMLl.MOD] 


#**»*«■**«■•«■ 


a- «■ 

* AR * 

1 *■=- 

«• *=s= 
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S2 MC =IX 
I 

I 

II 12 
I 

I VB AR 

I 
1 
I 

I 

II 22 
I 

I 

I 

I 

I 

I 

II 32 
I 

I 

I 

IL PT =S4 


* 


* IH HA 

«• 

MC 

«•< : 

* 

41 

* 

* 


«■ 




I 

I 

I 

I 

I 

113 

I 

=VB I 
V 


* «• 

* HA * 

«• 23 * 

«■ * 


**■***#**#«• 

I 

1 

I 

I 

133 

I 

1 

I 

I 

I 

43 


VB AR =VL «• # 

* PT » 

«• * 

I 

I 

I 

1 

14 115 

I 
I 
I 
I 
I 
I 

24 125 

I 
I 
I 
I 
I 
I 

34 135 

I 
I 
I 
I 

I 

44 4 5 


53 


52 
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Fig. 5.9 Schematic Diagram of [SIMLl.MOD] 
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TITLE=IMASA(LARGE GIGNAL) 

*TITLE = SOLAR ARRAY SYSTEM WITH SHUNT REGULATOR * 

* AND CONSTANT POWER LOAD « 

PARAMETER VALUES 

R AR = .001. C AR = 5E-4. L AR = lOE-6. LL AR = 1. 

LLSAR = 0. . TA AR = 301. , FC AR = SE-5, FV AR = -2E-3 

NP AR = 324. NS AR = 58 

■«••«••»■•»«■#•«••»*•»•«■«••«• «•«■«•*•« ■«•■«■*■«•#«• 

INITIAL CONDITIONS 
VB AR = 28. 141 

IL AR = 29. 694 

* ii- * 

INITIAL CONDITIONS 
XI HA = 10. 678 

S2 HA =- 18. 499 

X2 HA = -21. 650 

PARAMETER VAL UES 
VR HA = 12. a 

LI HA = i.5, K1 HA = 0.4548685 

A1 HA = -6. 00652 

Z1 HA = 2280. . PI HA = 2900. 

A2 HA = 1. , 72 HA = 667. . P2 HA 2900. 

20 HA = 539. 1998. PO HA == O. 

VF HA = 1 . 4 

K2 HA = 0.0150376, K3 HA = 0.083333 

L3 HA - 1.4E-3. Cl HA = 5450. , C2 HA -9.72 

I NAHA -- 40. 

*•»■**•«•*•«•* 

VR PT=15 

PWOPT=^ 0. SW PT^ 50000 
PW1PT=800 
TC PT=^5E-1.PC PT= j 
Cl MC=1,C2 MC= -1,C3 MC=1,C4 MC=0 
«• 

PRINTER PLOTS 

ONLINE PLOTS 

TMAX = 10E-3. TINC-=lE-6 

PRATE=100. OUTRATE=^ 50 

SIMULATE 

X I C — X 

PARAMETER VALUES 
PW2PT=^750. SW PT=^ 50000 
PW1PT=800 
TC PT=2E-3. PC PT=J 
Cl MC=1.C2 MC= --1.C3 MC=1.C4 MC=0 
«■ 

PRINTER PLOTS 
ONLINE PLOTS 
*INT M0DE=4 
DISPLAY 1 (OVERPLOT ) 

VB AR 

DISPLAY2 (0VERPL07) 

S2 MC. IH HA. IL PT 
TMAX==10E-3. TINC=lE-6 
PRATE=100. 0UTRATE=^50 
SIMULATE 


Fig. 5. 11 User's Parameter Values and Analysis Commands [ SIMLl . ANC ] 


119 


The simulation results are identical to those of the previous 
simulation, shown in Fig. 5. 7. 

For simulation of the Type-2 shunt regulator, the shunt 
regulator model [FS] replaces the Type-1 shunt model [HA] in 
the previous system configuration. This system model 
[SIML3.M0D] and schematic diagram are shown in Figs . 5 . 13-14. 
In this simulation (Fig. 5. 12), the constant power load de- 
creases linearly from 800W to 500W at a given rate (50W/ms). 
As the load current decreases, the shunt current increases 
accordingly. As a result, the solar array output current, 
again, becomes almost constant so the bus voltage remains 
within the specified voltage range. 

The difference between the bus voltage at 800W and the bus 
voltage at 500W results from the shunt current being a linear 
function of the bus voltage. The input data requirements and 
the analysis program (SIML3.ANC] used in this simulation are 
given in Figs . 5 . 15-16 . 
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V(BUS) CURRENT 

, 28.10 28.20 28.30 0.00 8.00 16.00 24.00 


8 



TIME 



Fig. 5. 12 Simulation Results of Type-2 Shunt Regulator Model 
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MACRO FILE NAME = MACROS 


* SOLAR ARRAY SYSTEM WITH TYPE 2 SHUNT * 

* REGULATOR AND CONSTANT POWER LOAD * 


LIST MACRO COMPONENTS = 
MODEL DESCRIPTION 
LOCATION = 1, AR, INPUTS 
LOCATION = 22, FS, INPUTS 
LOCATION = 6. PT, INPUTS 
LOCATION = 11. MC, INPUTS 
END OF MODEL 
PRINT 


AR. FS, PT 

= MC(S. 2=1 X) 

= AR(VB=VB) 

= AR(UB=VL) 

= FS(IH=S. 1 ). PT( IL=S. 4 ) 


Fig. 5. 13 System Model Generation Using Type-2 Shunt [SIML3.MOD] 


* AR * 

*=.-s:=s===:==s= = = = = =:==s==:=> 

A I 

S2 MC =IX I 
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Fig. 5. 14 Schematic Diagram of [SIML3.MOD] 
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TITLE=NASA(LARGE SIGNAL) 


* SOLAR ARRAY SYSTEM WITH TYPE 2 SHUNT REGULATOR * 

* AND CONSTANT POWER LOAD * 


PARAMETER VALUES 

R AR = .OOli C AR = 5E-4. L AR =-' lOE-6, LL AR = 1. 
LLSAR = 0. , TA AR = 301. , FC AR =- SE~5, FV AR = -2E-3 
NP AR = 324. NS AR = 58 
■»■■«•■«•■«■•«••»•«•«■»•■«•*«««•«■*«•«■*«•*■«■«■* 

INITIAL CONDITIONS 
VB AR = 28. 14 
IL AR = 30. 133 


■ft* *■«■*•«•■»** •«• * 

PARAMETER VALUES 

VR FS = 28.02. VSAFS = 15 

R1 FS = 12400, R2 FS =2. 2E6, R3 FS = 1.78419 
GM FS =4. VTHFS = 3, ILMFS = 45 

#«■#*■«■**■»• 

VR PT=15 

PW0PT==700. SW PT-50000 
PW1PT=800. PW2PT= 500. PW6PT = 500 
TC PT=5E-1.PC PT=^02 
Cl MC=1,C2 MCi=-l.C3 MC=1,C4 MC=0 
# 


PRINTER PLOTS 
ONLINE PLOTS 
*INT M0DE=4 
DISPLAYl 
VB AR 

DISPLAY2 (OVERPLOl) 
S2 MC. IH FS. IL PT 
DISPLAYS 
PW PT 

TMAX=10E-3, TINC=lE-6 
PRATE=tO. 0UTRATE=50 
SIMULATE 


Fig. 5. 16 User's Parameter Values and Analysis Commands 


[SIML3 . ANC] 


124 


5.2 PARTIAL SHUNT REGULATOR/SOLAR ARRAY MODEL 
5.2.1 INTRODUCTION 

The full shunt regulator has been widely used in space 
power systems. One disadvantage is that under a zero load 
condition it must dissipate the entire array output power. 
This complicates the thermal design and increases its size 
and weight. The partial shunt regulator reduces the peak 
power dissipation to 30 to 50 percent of the amount dissi- 
pated in a functionally equivalent full shunt. 

As in a full shunt, there is a control circuit that per- 
forms the same basic function and there are sections intended 
to dissipate or bypass some of the excess power. In a par- 
tial shunt system, the power sections are placed across only 
part of the total number of series-connected solar cells. 
Thus the shunt circuit is used to absorb part of the output 
of the solar array. The design of a partial shunt is greatly 
dependent on the type of solar cells used, the degradation 
encountered, and the temperature extremes to which the cells 
are subjected in orbital flight. 

As shown in Fig. 5. 17, the total number of series solar 
cells is expressed as 

Ng = N„ t Ng (5.7) 

cells 
cells 


where Ny = the momber of upper series 
N- = the number of lower series 
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ig.5.17 Partial Shunt Configuration 









The number of series cells, N , is selected to deliver a 

s' 

given fixed array's output voltage, = V^, throughout the 

life of the system [8]. 


N 


S 


d -V 
V mp 


V 


(5.8) 


where d^ is a voltage degradation factor equal to the ratio 
of the end-of-life to beginning-of-life voltage at the ar- 
ray's maximum power point, and at a steady-state array tem- 
perature. However, the array output voltage is lower than 
at the beginning-of-life if is chosen using Equation 
(5.8). This is not a proper region of operation since the 
system can easily become unstable. Therefore, the factor, a 
< 1, must be applied to to ensure the system is operating 
in the stable region. When the solar array output voltage 

for normal operation is greater than V , 

mp 


d • V = a -V 
S V mp r 


(5.9) 


In partial shunt configuration, the array is electrically 
split into two sections: upper and lower. The partial shunt 
element is placed across a part of the total number of series 
cells defined by the fraction m: 


m = + Ny ) 


(5.10) 
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When the shunt transistor is saturated ( = 0 assumed ), 

the highest possible voltage which the upper bank of Ny cells 

can deliver must not exceed V . This will occur when the 

r 

array is completely unloaded and the solar cells are at their 
lowest temperature, at the beginning-of-life . 

N..-V (max) < V (5.11) 

U oc'^ ' r 

By substituting Equations (5.7), (5.9) and (5.10) into 
equation (5.11), 

m > l-[o-d»V /V (max) ] (5.12) 

V mp ' oc' ' ' 

Thus, the location of the array tap-off point is seen to be 

dependent on the values of o, d , V and V (max). 

V mp oc 

5.2.2 MODEL GENERATION [AP] 

For the partial shunt configuration, as shown in Fig. 5. 18, 
another solar array model has been developed to provide the 
tap-off point of the solar array for the shunt regulator. 
The control voltage, V^, is generated by sensing the main bus 
voltage to control the shunt current. In this model, two 
solar array models are connected in series. The lower solar 
array supplies the shunt current to regulate the bus voltage. 
As a result, Newton iteration routines are included in the 
model and the output current of the upper solar array is 
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STATE EQUATIONS 
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Fig. 5. 18 Circuit Model of Solar Array/Partial Shunt (AP) 
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equal to the inductor current, The output current of 
the lower solar array is the sum of the shunt current and the 
inductor current. With the known upper and lower solar array 
output voltages, and V^ 2 ' ^wo state equations exist: 


= V 

'1 dt ^'lI 
dV„ 

' — — = I - I 

'1 dt ^L1 L 


where 


5.2.3 EXAMPLE OF SIMULATION 

simulation of the partial shunt model [AP], system 
IPSS5.M0D] (Fig. 5.20) is configured as shown in Fig. 5. 21. 
The system consists of the partial shunt model [AP], the full 
shunt model [FS] and a constant power model [PTj. As shown 
in the schematic diagram (Fig. 5. 21), the solar array model 
(partial shunt model) [AP] is connected with a full shunt 
model [FS] and a constant power model [PT] . 

As shown in Fig. 5. 19, as the load power decreases at a 
certain rate, the solar array output voltage increases. The 
bus voltage increases accordingly and becomes greater than 
(28.14V) which causes the shunt current to flow. As shown 
in Fig. 5. 18, the upper solar array carries only the load 
current and the lower solar array carries both the load cur- 
rent and shunt current. 
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V ( BUS ) CURRENT 
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TIME ixlO ms] 


Fig. 5. 19 Simulation Results of Partial Shunt System 
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* * ♦ * 


MACRO FILE NAME 


MACROS 


* PARTIAL SHUNT SYSTEM WITH CONSTANT POWER * 

* LOAD * 

■K- REVISED ON 3/3/1986 

LIST MACRO COMPONENTS = AP, FS. PT 
MODEL DESCRIPTION 

LOCATION--=li AP, INPUTS rS( IH=lH>, PT(IL=IL> 

LQCATI0N=--33. FS, INPUTS = AP(VB=VB) 

L0CATI0N=5, PT, INPUTS = AP(VB=VL> 

END OF MODEL 
PRINT 


Fig. 5.20 System Model Generation Using Partial Shunt [PSS5.MOD] 
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Fig. 5. 21 Schematic Diagram of [PSS5.MOD] 
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Since the load current continues to decrease, the upper solar 
array output voltage increases. The lower array current, the 
sum of the load and shunt currents, increases so the lower 
solar array output voltage decreases. The sum of lower and 
upper solar array voltages becomes the bus voltage. 

As a result, the bus voltage could be regulated within the 
specified range. The input data requirement list and analy- 
sis program used for this simulation are given in 
Figs . 5 . 22-23 . 
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INPUT DATA REQUIREMENTS LIST 
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Fig. 5. ’22 Input Data Requirements List 


TITLE=NASA(LARGF SIGNAL) 

«^TITLE = PARTIAL SHUNT SYSTEM WITH CONSTANT * 

* POWER LOAD * 

*#*«• •«•«•«■ «•«•***#**■«••«•*<{•*«■■«■**<{■*•«■■»•■«•■«•«■«••»!•«■***«• 

PARAMETER VALUES 

R AP = .001, C AP = 500E-6. L AP = IE -6 
NP AP = 324, NSIAP - 29, NS2AP =29, LL AP =1. 
ILMAP = 45, VIHAP = 15.976, VU AP = 12.57 
V2HAP = 15.976, V2LAP = 12.57 

INITIAL CONDmONS 
VB AP = 28. 161 
ILIAP = 28. 409 
•»**»*»** 

PARAMETER VALUES 

VR FS = 28. 14, VSAFS =15. 

R1 FS = 12400, R2 FS = 2. 2E6, R3 FS = 1.7B419 
GM FS = 4. , VTHFS = 3. , ILMFS = 45. 

VR PT=15 
SW PT= 150000 
PW1PT=800, PW2PT=0 
TC PT=5E-3, PC PT=2 

# 

PRINTER PLOTS 
ONLINE PLOTS 
*INT M0DE=4 
DISPLAYl (OVERPLOT) 

VB AP 

DISPLAY2 (OVERPLOT) 

ILIAP, 12 AP, IH FS 
DISPLAYS (OVERPLOT) 

VOIAP, V02AP 
DISPLAY4 (OVERPLOT) 

PW PT 

TMAX=10E-3, TINC=lE-6 
PRATE=200, 0UTRATE=50 
SIMULATE 


Fig. 5. 23 User's Parameter Values and Analysis Commands [PSS5.ANC] 
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CHAPTER 6 


BATTERY-DISCHARGER MODELING 


6.1 INTRODUCTION 

The basic functional requirement of the battery discharger 
is to regulate the bus voltage when the solar array is inca- 
pable of meeting the load demand. In this mode both the solar 
array and the battery function as power sources. Figure 6.1 
shows the system operating in the battery discharge mode. 
When the bus voltage drops below a certain magnitude, the 
control signal activates the battery-discharge converter in 
order to maintain regulation of the bus voltage. 

In the Direct Energy Transfer (DET) system, when the solar 
array output power decreases due to the weakening of the il- 
lumination level, the battery discharger operates from no- 
ibad condition through the converter 

discontinuous-conduction mode (DCM) of operation and ulti- 
mately to the continuous-conduction mode (CCM). Since the 
steady- state and dynamic behaviors of switching regulators 
for these two modes of operation are quite different, the 
design of the feedback controller must take 
atioh the two modes of operation. 
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into consider- 




6.2 BATTERY-DISCHARGER MODELING AND ANALYSIS 

The power stage of the battery-discharge-converter used 
in the DET system is shown in Fig. 6. 2. The circuit basically 
operates as a buck-type regulator operating in a push-pull 
mode. Since the terminal voltage of the battery is between 
17 and 25 Volts, an autb-transformer is used to boost the 
voltage to 28V. In this study, a simplified circuit is used 
as shown in Fig. 6. 3(a). The two sources in Fig. 6. 3(a) rep- 
resent the transformer action between the power switches' ON 
and OFF intervals. 


6.2.1 ANALYSIS OF BATTERY- DISCHARGE CONVERTER IN DCM 

As illustrated in Fig. 3. 2, there exist three intervals, 
d^T , d«T_ and d_T . While the first interval, d.,T , 
(inductor charging) is dictated by the "ON" time of the 
switching function and is a known quantity, the second in- 
terval, d>,T , (inductor discharging) is unknown and depends 
in general bn both the length of the first interval and cir- 
cuit parameters. 

Using ideal elements in the circuit in Fig. 6. 3 (a) and ig- 

nbring load dynamics, three switched networks in the 

discontihuous-cbnductibn mode are shown in Fig. 6. 3(b). For 

T 

the choice of the state vector, x = [i v] , the state-space 
equations bf the three switched networks become: 


13^8 


Output Current 
sense 
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‘Txagure 6.2 Battery discharger converter circuit 





X 


for interval d,T 
1 s 


= AjX + 

X = A„x + B,v for interval d_T (6-1) 

X = A^x + B-v for interval d_T 

o 3 g 3 s 

where 



( 6 . 2 ) 


By averaging equation (6.1), the following linear, alge- 
braic, static equations result. 


0 1 


I 


D^(1+k)+kD2 

L 



+ 

L 

- Jl 

C CR 


V 


0 


A 


X 


B 


Then the steady-state solution of Eq. (6.3) is 


~ . M = 

g 


D.,(l+ic) + kD, 


I = 


R (D^ + D^) 


(6.4) 


Hence, the dc conditions depend not only on and R, as in 
the continuous-conduction mode, but also on D 2 • Using the 
additional constraint on the average of the inductor current 
for (D^^D 2 )T 3 , 


I 


(l+K)Vg - V 
2L 


D^T 
1 s 


(6.5) 


The dc conditions are determined with the known quantities 
and R. For example, substitution of Eq.(6.5) into 
Eq.(6.3) results in 


D 


2 


1 

T 



+ [( 



+ 4t^(1+k)]^ 


6 . 6 ) 


where 

T = 2L / (R T ) ; normalized time constant. 
XI s 
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The boundary condition between the two conduction modes is 


D 2 < 1 - ; discontinuous conduction mode 

D 2 > 1 - ; continuous conduction mode (6.7) 

Using Eqs.(6.6), the critical value of can be determined 
by solving D 2 =l-D^. 

^ ■ ^1 

^ncrit 1 + K + k(1 - D^)/D^ (6.8) 


Suppose that the output voltage is regulated with a proper 

feedback circuit, then the dc voltage gain, M, can be treated 

as a constant for a fixed input voltage. The duty ratio, D^, 

in the discontinuous-conduction mode can be derived as a 

function of M, or as a function of t . 

n 


^1 = 


T ( M- K ) 

1 + K - M 


(6.9) 


Using the results in Eqs.(6.4) through (6.9), the dc be- 
havior of the battery discharger is illustrated in Fig. 6.4. 
In Fig. 6. 4(a), the boundary condition, '^jj^rit' 
two modes is plotted as a function of The dc voltage 

gain, M, as a function of for different values of x^, is 
shown in Fig. 6. 4(b). 
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6.2.2 COMPENSATOR DESIGN FOR BATTERY DISCHARGER 


Since the battery discharger operates from no-load to 
full-load, it is important to understand the start-up behav- 
ior of the converter. Suppose now that the solar array out- 
put power gradually decreases. Initially the 
battery-discharging current is zero (i.e., '^^^=0) and the duty 
ratio D^=0. As the solar array output power decreases, the 
operating point moves along the solar array curve which 
causes the bus voltage to drop. This activates the battery- 
discharger electronics via an error amplifier. In this 
start-up mode, it is very important to design a proper error 
amplifier circuit for the following reasons. 

Since the sensitivity of the dc voltage gain, with respect 
to the duty ratio, is very high in the light-load condition 
(start-up mode), as shown in Fig. 6. 4(b), any delay in the 
.logic circuitry or the false triggering due to noise will 
cause a high-frequency oscillation in the bus voltage. While 
in a low-gain feedback system the controller may not act 
quickly enough to regulate the bus voltage, since the plant 
time constant of the output capacitor and load resistance is 
very slow. This may cause a large-amplitude, low-frequency 
oscillation. This is especially critical when the output 
voltage of the battery-discharge converter is set such that 
it operates close to the solar array maximum power point. 
Since the low-gain controller allows the bus voltage to drop 
before the controller starts correcting the error, it is 
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possible that the system may move into the unstable region 
as discussed previously. Since the battery discharger oper- 
ates from the discontinuous-conduction mode to continuous- 
conduction mode, the error amplifier should be able to 
stabilize the system in both modes. 

Small-signal analysis of the battery discharger in 

discontinuous-conduction mode shows that the transfer func- 

fgl 

tions have a single pole and no zero . Any converter that 
can be described with two state variables, one of them being 
an inductor current, has a single-pole transfer function in 
the discontinuous conduction mode. Thus, the dynamic behav- 
ior of the converter, operating in the discontinuous- 
conduction mode, is completely different from that of the 
continuous-conduction mode where the system has a complex 
pole pair. Therefore, the feedback controller should com- 
pensate the plant dynamics in both modes of operation. Fig- 
ure 6.5 illustrates the root-locus diagram of the open-loop 
converter as a function of load resistance. Without a 
compensator, the system eigenvalue is close to the origin, 
thus, the system response should be very slow in the 
discontinuous-conduction mode as predicted in the dc analy- 
sis. Decreasing the load resistance beyond 'tbe con- 

verter will operate in the continuous-conduction mode, and 
the open- loop eigenvalues become complex. In this range, the 
response is slow and oscillatory. A further decrease of the 
load resistance will move the locus toward the circle with 
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Figure 6.5 Root locus of the open loop buck-type converter as 
a function of load resistance, R for DCM and CCM 
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the radius of = 1//LC, and finally approaches to the real 
axis again. 

Figure 6.6 illustrates the root-locus diagrams for dif- 
ferent compensators in both DCM and CCM. Let us first con- 
sider the converter operating in DCM. Since the open- loop 
characteristics have a single pole near the origin, to 
speed-up the response a simple proportional gain controller 
can be used to move the eigenvalue away from the origin. This 
is shown in Fig. 6. 6(a) . However, this gain also moves the 
eigenvalues to the imaginary- axis in the CCM operation. A 
lower gain may stabilize the system in CCM, but the slow re- 
sponse in DCM may cause a severe, bus voltage droop during 
transient. 

A typical integrator and lead-lag compensator as illus- 
trated in Fig. 6. 6(b) is examined next. The integrator is 
.used for the tight regulation of the output voltage. The zero 
compensates the excessive phase lag due to the presence of 
an integrator, and the pole attenuates the high-frequency 
switching noise. Since the integrator slows down the system 
response in DCM, the roots stay near the imaginary axis in 
DCM. 

Therefore, to achieve an optimum response for both CCM and 
DCM operations, a compensator must include a proportional 
gain to speed up the system response in DCM and the lead- lag 
type compensator for stable operation in CCM. This is il- 
lustrated in Fig.6.6(c). 
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Figure 6.6 Root loci of the closed loop buck-type converter 
in DCM and CCM for three compensators 

(a) proportional gain 

(b) integrator plus lead- lag 
Cc] proportional plus lead -lag 


compensation pole 
• compensation zero 


IL 



6.3 EASY5 MODELING FOR THE BATTERY- DISCHARGE CONVERTER [BP] 

In the DET system, as shown in Fig. 6.1, the output 
capacitor of the battery discharger is shared with the 
output-filter capacitor of the solar array. Thus, the 
capacitor is not included in the power stage model. Since 
only one state variable, inductor current is needed in 
the power- stage model, the state equation for each interval 
is as follows. 


For D^T (switch is on). 
Is' ' 


ij^ = [(1 + K)v^- i^R^- V^l/L 


For D«T (switch is off, i- 0), 
z s u 


= ■''r 


For D_T (switch is off, i_= 0), 

O S Li 


° 


For 

Chapter 


the feedback control, 
3 can be used. 


macro modules developed in 
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6.4 BATTERY DISCHARGER MODEL SIMULATION 

To simulate the performances of the battery discharger, a 
simplified DET system is configured using the following mod- 
ules . 

• [AR] solar array 

• [FS] full shunt regulator 

• [BP] battery-discharger power stage 

• [PT] constant power load 

• [WM] PWM( constant frequency control) 

• [MC] summer{EASY5 standard module) 

• various compensator modules [MP] , [PZ] , [DG] 

Figure 6.7 shows the EASY5 generated system configuration 
[SIM.MGLj. The following simulations are done with the 
single-loop control employing constant- frequency PWM. Two 
compensators are employed, as shown in Fig. 6. 8, to verify the 
previous analysis. Figure 6.9 shows the step- load transient 
response from 1300W to 1400W when the battery discharger op- 
erates in the continuous-conduction mode. As shown in the 
figure, the system behaves equally well for both 
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Figure 6,7 EASY5 generated system configuration 
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compensators. In Fig. 6. 10, the load power is gradu- 
ally increasing and various modes are illustrated. The sys- 
tem initially operates in the shunt mode until 0.6 ms, then 
it gets into the dead-hand mode where the solar array solely 
supplies the load power. The operating point moves along the 
solar array curve as the load power increases. The battery 
discharger, employing the proportional-gain feedback 
(Fig.6.6(c) ) , starts regulating the bus voltage at t=lms as 
shown in Fig. 6. 10(a). While using the integrator feedback, 
it starts to operate at t=2ms and slowly creeps up to reach 
the output-voltage regulation, as shown in Fig. 6. 10(b). As 
shown in this figure, the slow response in DCM of the 
integrator- type controller causes severe droop in the bus 
voltage. 

The results in Figs. 6. 9 and 6.10 show that a compensator 
designed for the CCM operation may not be suitable for the 
DCM operation, and consequently, the compensator design for 
the boost converter must be done more carefully. This is 
because such a converter has a single pole in the 
discontinuous-conduction mode, and a complex-pole pair and a 
right-half-plane zero in the continuous-conduction mode. 
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CHAPTER 7 


SYSTEM LEVEL MODELING, ANALYSIS AND SIMULATION 
7.1 INTRODUCTION 

Solar arrays, as the primary power source for spacecraft, 
generate electrical power proportional to the illumination 
level. The operating point of the solar arrays is determined 
by the power conditioning equipment and the load character- 
istic. The power conditioning equipment consists of a shunt 
regulator, a battery bank and its charging and discharging 
electronics. When the output power of the solar array ex- 
ceeds the load demand, a portion of the excess power is used 
to charge the battery bank, the secondary power source, and 
the rest dissipated in the shunt regulator. When the load- 
power demand is higher than the solar array output power, the 
battery bank supplies the additional load-power demand 
through the battery discharging electronics. Since the solar 
array output power varies over a wide range according to the 
changes in the illumination level, several modes of operation 
exist depending on the characteristics of the solar array, 
the power-conditioning equipment, and the load. 

As described in Chapter 4, the solar array output charac- 
teristic is nonlinear. This nonlinear-source coupled with 
nonlinear load characteristics may result in multiple equi- 
librium points under a given operating condition. The actual 
equilibrium state that the system resides on is determined 
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by the stability nature of the equilibrium points and its 
past history. Thus, it is necessary to analyze the stability 
of the equilibrium points and to understand various modes of 
operation. 

In Section 7.2, stability of the equilibrium points of a 
solar array system with a nonlinear load is analyzed. The 
behavior of different solar array operating points are plot- 
ted as trajectories on the state plane. Stability and tran- 
sient responses of the system operating near the solar array 
maximum-power point is considered in Section 7.3. 

In Section 7.4, various modes of operation of the DET 
system are described. The solar array switching system and 
the partial shunt system with a buck converter are simulated 
in Section 7.5. Conclusions are presented in Section 7.6. 
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7.2 LARGE- SIGNAL BEHAVIOR OF THE SOLAR ARRAY OPERATING POINT 

Fig. 7.1 illustrates a solar array power system, in which 
the static load line represents a typical regulator load- 
line characteristic. The voltage Vq in the figure is the 
minimum solar array output voltage above which the converter 
is able to regulate its output voltage. When the converter 
is regulating, the solar array sees a constant power. Note 
that the converter exhibits a negative-resistance character- 
istic in this mode. Below the converter behaves as a 

resistive load to the solar array. 

As shown in Fig. 7.1, there exist three equilibrium points 
in the system. The large-signal behavior of a solar array 
operating point depends on the stability nature of each 
equilibrium point. The qualitative large-signal analysis can 
be carried out by considering a second-order system as shown 
in Fig. 7. 2. The solar array output and the load character- 
istics are replaced by nonlinear functions. The LC low-pass 
filter in the model represents actual physical elements in 
the system such as an input filter of the switching regulator 
or the cable impedance between the solar array and the load. 

Let us first consider the stability nature of each equi- 

T 

librium point. Defining the states x = [ij^ system 

state equations are 


dt 



1 

— ( -f(iL) - ) (7.1) 

Li 
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Figure 7.2 Equivalent circuit model of the system in Fig.l .1 
including LC fitter 
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1 


1 


dt 


= “ ( - iR ) = - < > ("-2) 


where v<, and are defined in Fig. 7.2. 


From Eqs.(7.1) and (7.2), the equilibrium points in Fig. 7.2 


are 


(7.3) 

II = f(V^) (7.4) 


The behavior of the trajectories in the neighborhood of an 
equilibrium point can be found by linearizing Eqs.(7.1) and 
(7.2) about an equilibrium point. The nature of these equi- 
librium points is analyzed by introducing small perturbations 
about the equilibrium points. 





(7.5) 


V 


C 


" ^c 


(7.6) 


where and represent the steady-state quantities, and 
ijj^ and are the small- signal perturbations. 

Substituting Eqs. (7.5) and (7.6) into (7.1) and (7.2) results 
in 



dt 


T ‘ - < ''c * ''c > ' <’■’> 

Li 
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1 


(7.8) 



dt C 
Expanding Eqs . (7.7) 
equilibrium points, 
imation are 


[ -f(V(, ♦ v^) - ( ^ ) 1 

and (7.8) in a Taylor series about the 
the equations of the first order approx- 




(7.9) 


where 




df(ij^) 


di. 


and similarly. 


f(V^ + v^) = f(V^) 



(7.10) 


where 


1 

r, “ 


dv. 


where r., and r^ are the incremental resistances (tangential 
slope) of the solar array output curve and the load line, 
respectively, at the equilibrium points. 

Employing Eqs. (7.3) through (7.10), Eqs. (7.7) and (7.8) can 
be simplified as 

di^ ( 1/L ) ( -rg i, - ) 

dv^ ( 1/C ) ( i^ - (1/rj.) 
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Note that the value of Tg is the same as the output resistance 
at a point on the solar array I-V curve. The system 
eigenvalues are expressed as 

1 1 4 1 

-( + ) + [( + )" - (1 + ) 1 ^ 

L rj^C L r^C LC r^^ 

(7.12) 

From Eq.(7.12), the stability of each equilibrium point can 
be determined for given parameter values. For the load line 
shown in Fig. 7.1, the following results are observed. 

In Region ttl : rj^ > 0 

The equilibrium point, is stable and is either 

"node” (real) or "focus" (complex). 

In Region #2 : r^^ < 0, |rg| < |rj^| 

and X 2 are both real but of opposite sign, and Point 

P is a "saddle" and is unstable. 
e2 

In Region #3 ; r^^ < 0, | rg | < |rj^| 

The stability of Point P^g depends on the parameter 
values . 

Since the normal mode of operation should be in Region #3, 
the values of L and C in Fig. 7.2 must be chosen to provide 
the stable equilibrium points in this region and to allow 
further analysis. However, it must be noted that the equi- 



163 


librium points in Region #1 are always stable and those in 
region #2 are always unstable (saddle points) regardless of 
the parameter values. The value of rj^ in the constant power 
region is 




3i 




ir=lR (7.13) 


Approximate values of r^, and r, can be found simply by draw- 

O 1-1 

ing tangential slopes to the solar array I-V curve and load 
line, respectively, for each equilibrium point. 

To construct a state-plane trajectory the isoclines for 
each equilibrium point are derived using Eq.(7.11). 


di 


£ 


dv 

c 


= 0 



(7.14) 




dv 


c 



(7.15) 


The results in Eqs.(7.14) and (7.15) are the tangential 
slopes to the source and load curve at each equilibrium 
point, thus the trajectories cross the solar array curve 
horizontally and the load line vertically. Near the equi- 
librium points which have real eigenvalues, the trajectories 
follow the slope of the eigenvectors of the system, Z. 
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^2 = _ 


r 


SI "SI 4 1 

_) + [(—- )" 1 '^ 

(7 

L r^C L r^C LC 


16) 


Using the above results, the state plane trajectories can be 
sketched. Fig. 7. 3 illustrates the trajectories for the case 
where the equilibrium point is a stable node. As shown 
in the figure, the trajectories are separated by the 
separatrix that has the slope at Pg 2 • starting 
points and ^2 7.3 approach two entirely different 
steady states, even though the initial states differ by only 
a small amount. 

For the solar array model used in this study the following 
cases are set up for the computer simulations. As mentioned 
above, if we assume that point P^^ Region #3 is stable then 
there can be four possible cases: 


i) 

^1 

- node. 

^e3 

- node 

ii) 

^1 

- node , 

^e3 

- focus 

iii) 

^el 

- focus. 

e3 

- node 

iv) 

^el 

- focus. 

^ -5 

e3 

- focus 


However, investigation of the cases i) and ii) are sufficient 
to understand the behaviors of the system. The results from 
simulating these cases are shown in Fig. 7. 4(a) and 7.4(b) 
given the following system operating conditions. 
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ITRIX 
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The solar array short-circuit current, ^gg/ is 45.41A, the 

open-circuit voltage, is 32.0V, and the maximum power 

is 1050 watts. The constant-power load line is set to 850 

watts, and the boundary of Region #1 and #2, V„, is 15V. The 

K 

values of rg and r^^ at the equilibrium points in each region 
are listed in the following table. 



^el 

^e 2 

P 0 
e3 


45.38 

45.02 

30.34 


12.03 

17.31 

28.02 


64.92 

37.3 

0.20 


0.025 

-0.56 

-1.47 


The cases i) (^^3 node) and ii) (Pq 3 focus) are shown in 
Fig. 7. 4a and 7.4b. The eigenvalues and the slope of the 
eigenvectors for the particular values of L and C are listed 
in the figure. In Fig. 7. 4, the separatrix is clearly illus- 
trated. It is obtained by adjusting the initial condition 
that is located near the line estimated from information 

The trajectories from the initial points 
in the left-hand side of the separatrix (P^ through Pg) con- 
verge to Pi/ and those in the right-hand side of the 
separatrix (P.^ through P^ 2 ) ^e 3 ’ 


about ^2 


Since in Region #1 the system bus voltage is not regu- 
lated, the desirable mode of operation should be in Region 
#3. This implies that there needs a limiting mechanism that 
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Figure 7.4(a) Simulation of state- plane trajectories for the 

case where the equilibrium point is 'stable node' 

(L = 10 ]iH, C = 500 vF) 
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Figure 7.4(b) Simulation of state-plane trajectories for the 

case where the equilibrium point is 'stable focus' 

[L = 30 vH, C = 500 vFj 
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prevents the operating point from settling into Region #1 or 
a triggering method to move it to Region #3 even if the op- 
erating point starts in Region #1. The limiting mechanism 
will be described in the following section. A triggering 
method is illustrated as follows. 

Suppose that the operating point is at rest in Region #1. 
One way to move this operating point into Region #3 is to 
decrease the load power below the solar array curve so that 
the only one desired equilibrium point exists. For the load 
line shown in Fig. 7.5, the equilibrium points associated 
with Region #1 and #2 have moved outside these regions, and 
therefore they no longer exist. However, the behavior of the 
state trajectories can be still seen by defining virtual 
equilibrium points, and which the slopes of 
the eigenvectors are also sketched. 

Fig. 7. 6 shows the simulation result using the previous 
system for the case i). In Fig. 7. 6(a), the trajectory starts 
from an arbitrary initial point, and moves to 1^2 for a load 
power of 800 watts (Curve #1). As the load power contin- 
uously decreases, the operating point moves along the solar 
array curve until it reaches Curve #2. Further decrease in 
the load power beyond this point moves the operating point 
from to P^, then to P^ for Curve #3. Fig. 7. 6(b) shows the 
time history of v corresponding the trajectory from P., to 

C J. 


170 



Figure 7.5 A trajectory from the initial operating point 
in Region §1 to Region #2 
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TIME ""see 


Figure 7 . 6 Simulation of the triggering method from the 
undesired stable region to the desired region 
fa) state- plane trajectory 
(b) time history of the bus voltage 
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The above stability analysis for the system in Fig. 7. 2 can 


be extended to the following cases. 


System in the battery discharging mode : Figure 7.7 il- 
lustrates the system in the battery discharging mode of op- 
eration. The battery discharging current can be expressed 
as a function of the bus voltage with an output admittance, 
l/r„. The state equations are then. 


dt 




(7.17) 


dv 


C 


dt 


1 

C 


iL - i - f(^c> ) (7.18) 

B 


Following the same procedures as in the previous derivations, 
the state equations near the equilibrium point are 

di, ( 1/L ) ( -rg i, - ) 

dv^ ( 1/C ) [ i^ - v^ ] 

(7.19) 
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Figure 7.7 Simplified block diagram of the system in Fig. 7.2 
with the shunt regulator or the battery discharger 


{ Y , : for shunt mode 

ifr„ : for battery discharge mode 

D 
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Then the eigenvalues are obtained as 



4 

LC 


(1 + 



(7.20) 


where 

■^E = I 

Comparing Egs.(7.12) and (7.20), the additional term, l/r„, 

x5 

stabilizes the system, which can be explained as the system 
is operating on the stiff voltage source line. 


System in the shunt mode ; Figure 7.7 also can be used 
to illustrate the system in the shunt mode of operation. In 
this mode of operation, Eqs.(7.17) and (7.18) are valid if 
the 1/^g term is replaced by the shunt transadmittance, 
Figure 7.8 shows the state plane trajectories and the time 
history of the states, with and without the shunt regulator. 
This figure clearly shows the effect of the additional term, 
l/r„ described above such that the system without the shunt 
regulator is quite oscillatory while that with the shunt 
regulator has faster responses and the system appears to have 
real valued eigenvalues. 
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7.3 ANALYSIS OF THE OPERATING POINT NEAR THE SOLAR ARRAY 
MAX I MUM- POWER POINT 

To utilize the maximum available power from the solar ar- 
ray, several approaches were proposed in [ 7 ] and [10] to 
place the operating point as close to the solar array 
maximum-power point as possible. However, near the solar 
array maximum power point, the dynamic resistance of the so- 
lar array is very sensitive to changes in the operating con- 
dition and, consequently, may result in system open loop 
instability. 

To derive the solar array maximum-power point, one can 
differentiate the constant-power load line, P = VI, and set 

cL 

it to zero. 

I = 0 (7.21) 

(7.22) 

Equation (7.22) states that at the maximum-power point, the 
constant-power load line, r^^, is equal to the dynamic output 
resistance of the solar array, r^, as shown in Fig. 7.9. 


3P 31 

± = V + 

3V 3V 

V 3V 

I 31 
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Figure 7 . 9 Relationships of solar array characteristic and constant 
power load line at the maximum power point, M 
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To analyze the behavior of an equilibrium point near the 

maximum power point, one can consider the solar array small- 
[ q 1 

signal model. ^ One of the system eigenvalues at the equi- 
librium point is 


X 


1 


: ( 1 + R / r^ 

eq ' eq ' L 


(7.23) 


where is the equivalent shunt capacitor in the solar ar- 

ray small-signal model. 


From Eq.(7.23) , for 


o 

V 

the 

system is: 

stable 

if 

1 ! 
' o ' 

1 < 

1 r^ 1 , and 

unstable 

if 

|r 
' o 

1 > 

I^l' 


In the stable region near the maximum-power point, M, the 
magnitudes of r^ and r^^ are close together and the eigenvalue 
is close to the origin, which implies very poor transient 
response and stability margin. In this respect it is desir- 
able to operate the system in a region where |r | << | r^ | . 
This can be seen from the step transient response simulation 
using the system model developed in Section 7.2 (Fig. 7.2). 
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Figure 10(a) shows the response of the step change in the 
load power from 800 watts to 790 watts, and Fig. 7. 10(b) from 
1040 watts to 1030 watts which is close to the maximum-power 
point (1050 watt). Therefore, to design a controller to 
utilize the maximum available power from the solar array, one 
must take into consideration of the stability of the operat- 
ing point and the transient behavior of the system. One way 
to overcome this problem is to use the shunt regulator to 
dissipate some solar array power. In this mode of operation, 
the system can be stabilized as shown in Section 7.2. 
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(D 



8 



Figure 1.10 Step load transient responses 

(o) 800 to 790 wott step 
(b) 10U0 to 1030 watt step 
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7.4 ANALYSIS AND SIMULATION OF THE SYSTEM OPERATING MODES 


A direct energy transfer (DET) system whereby the primary 
power source, a solar array, is coupled through a main dis- 
tribution bus directly to the spacecraft electrical loads. 
The various power-conditioning components, as shown in 
Fig. 7. 11, are activated only as needed, thus requiring the 
system to process only the amount of power needed to maintain 
the bus at the specified voltage level. A system mode of 
operation is determined by the power system control unit. 
The Central Control Unit (CCU) continuously monitors the bus 
voltage and commands the system to operate in the battery- 
discharging mode when is less than the battery 
discharger regulation voltage, or in the shunt-active mode 
when is greater than the shunt regulation voltage. 
When V, , is between V , and V , there exists a "dead band" 
at which the load power is equal to solar array output power. 
When the shunt is activated, any variation of the bus voltage 
generates a proportional shunt current, The shunt cur- 
rent together with the load current sinks the solar array 
current and regulates the bus voltage. When the system op- 
erates in the battery-discharge mode, the bus voltage is 
regulated by the battery discharger. The solar array current 
is set by the bus voltage on its output I-V curve and the 
battery-discharging current is the difference between the 
solar array current and the load current. 
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CABLE /I B load I CABLE /I ^bus 
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Figurel .11 Direct Energy Transfer [DET) spacecraft power system 





In this work, modeling of the battery and the battery- 
charging electronic are omitted, thus, the battery is treated 
as an ideal voltage source. The entire system model includes 
solar array, shunt regulator, battery discharger and 
constant-power load. In order to simulate the system mode 

of operation including transition mode between shunt regula- 
tor and battery-discharging modes, a simplified DET system 
is configured using following component modules. 

• [AR] — •- solar array 

• [FS] shunt regulator 

• [BP] battery discharger power stage 

• [PZ] two-pole two-zero compensator 

• [WM] PWM( constant frequency control) 

• [PT] constant power load 



Figure 7.12 [SIM. MOD] displays user input model description 
data. Figure 7.13 [SIM.MGL] shows EASY5 generated DET system 
configuration. Figure 7.14 [SIM.MGL] is EASY5 generated re- 
quired input data list, and user furnished input parameter 
values and simulation commands are shown in Fig. 7. 15 
[SIM.ANC] . 
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*•«•****•«• it * -R- * * 

DET(NASA) SYSTEM LARGE SIGNAL MODEL ituprs- 
REVISED 3/-1/86 

•»■ it -R- «*•«•»« #■)«••»•«• it «•#«•* it# it »•»■«■•«■#«■■«■«»■«•■»•■«•*«■«•■«• #*it it it # 

MACRO FILE NAME^-MACROS 

LIST MACRO COMPONENTS=AR, FS, PT, BP, PZ, WM 

MODEL DESCRIPTION 

LOCATION=4. Af<> JNPUTS=MC(S, 2=1 X) 

L0CATI0N=25, FS, INPUTS=AR ( VB=VB ) 

LOCATION=l, PT, 1NPUTS=AR(VB=VL) 

L0CATI0N=43, BP, 1NPUTS=AR ( VB=V2 ) , WM( I Q= IQ ) 

L0CATI0N=16, PZ, 1NPUTS=AR ( VB=VO) 

L0CATI0N=46, WM, ] NPUTS=PZ < VE=VE ) 

L0CATI0N=22, MC, INPUTS=FS( IH=S, 1 >, BP( ILL=S, 3), PT( IL=S, 4) 

END OF MODEL 

PRINT 


Figure 7,12 User Input Model Description 
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Figure 7 . 13 EASY5 generated system configuration 
[SIM.MGL] 
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TITLE=NASA<LARGE SIGNAL) 

PARAMETER VALUES 

R AR = .001, C AR = lOOOE-6, L AR = lE-6. LL AR 1. 

LLSAR = 0., tA AR = 301. > FC AR = 8E -5, FV AR = -2E-3 

NP AR = 324. NS AR = 58 

INITIAL CONDITIONS 
VB AR = 28. 141 
IL AR = 29. 694 

*»FULL SHUNT 

PARAMETER VALUES 

VR FS 28. 14. VSAFS = 15 

R1 FS = 12400. R2 FS =2. 2E6. R3 FS = 1.78419 
GM FS =^4. VTHF3 = 3. ILMFS = 45 

*#*■»•■«••«•»» 

PARAMETER VALUES 

Cl MC=1.C2 MC= -1>C3 MC=1,C4 MC=1 

PC PT=1 
VR PT-15 

PW0PT=700. SW PT=50000 
PW1PT=700. PW2PT^i000 
TC PT=1 

INITIAL CONDITIONS 
IL BP-O 

PARAMETER VALUES 
L BP=100E-6 
RL BP=lE -3 
VI BP=20 
V2LBP=28 
GAMBP-. 75 
EPSBP-lE-3 

K PZ=. 2153625. ER PZ=6 
*WM PZ=23. 3E3 
WM PZ=^43E3 

WZ1PZ = 1256. WZ2PZ=6. 28E3 
WP PZ=50E3 

TI WM=20E-6. VP WM=3. VQ WM=. 5 

ER WM=6, CICWM=0, SCMWM=0 

VCXWM=10. . VCNWM=0 

DMXWM=1, DMNWM=^0 

ILXWM=20 

PRINTER PLOTS 

ONLINE PLOTS 

I NT M0DE=4 

D I SPLAY 1 

VB AR 

DI SPLAY2 < OVERPLOT ) 

IL AR. IH FS. IL PT, ILLBP 
DISPLAYS ( OVERPLOT ) 

VR WM. VC WM 
TMAX=1. 2E-3, TINC= IE -6 
PRATE=100. 0UTRATE=10 
SIMULATE 
X I C— X 

PARAMETER VALUES 
«TC PT=. 2E-3 
LLSAR=--250 

TMAX=1. 99E-3. TINC=2E-7 
SIMULATE 


Figure 7.15 User Input Analysis Program [SIM.ANC] 
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The system operating condition is subject to change due 
to a constantly varying load- power consumption and the solar 
array power generations. The DET system shown in Fig. 7. 11 
is simulated to observe how the system operating point moves 
from one mode to another. Figure 7.16 illustrates the case 
in which the solar array I-V characteristic varies according 
to the illumination level. In this example, the illumination 
level is set to be linearly time varying from Curve #1 (shunt 
active mode) to Curve #4 (battery discharging mode) as il- 
lustrated in Fig. 7. 16. 

From Fig. 7. 16 and the simulation result in Fig. 7. 17, the 
system mode of operation is described as follows: Starting 
from Curve #1 and the corresponding operating point A, the 
operating point travels along the shunt regulation line, V 

iT S 

until the solar array power decreases to Curve #2. In this 
mode, the shunt regulator is activated. It absorbs the ex- 
cessive current generated by the solar array through the 
shunt elements Between Curve #2 and #3 the dead band 
mode exists where both the shunt regulator and battery 
discharger are deactivated. Thus, the operating point trav- 
els along the load line from point B to C. Further decrease 
in the illumination level activates the battery discharger. 
The operating point then follows the battery-discharging 
regulation line, In this mode, the load power is sup- 
plied by the solar array and the battery. 
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V 


Figure 7.16 DET system mode of operation for the solar array 
illumination level changes 
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Figure 7.18 illustrates the case where the load power is 
linearly varying while the solar array I-V curve is fixed. 
From the Fig. 7. 18 and the simulation results in Fig. 7. 19, the 
system mode of operation is described as follows: With the 
solar array curve fixed and the load line below Curve #2, the 
operating point is set at Point A. As the load power in- 
creases beyond Curve #2, the shunt regulator is deactivated, 
and the operating point travels on the solar array curve un- 
til Point B. Further increase in the load power beyond Curve 
#3 activates the battery discharger, and the operating point 
moves along the battery discharger regulation line to Point 
C. 

Figure 7.20 shows simulation results when load is changed 
abruptly at 0.2ms from curve #1 to curve #4 in Fig. 7. 18. As 
the solar array can not afford the demand of load power, the 
shunt regulator is deactivated promptly and battery takes 
over the shortage. 








Figure 7 . 18 DET system mode of operation for the load power changes 
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1.20 1.60 2.00 


0 ^ fi 
u e 



0.00 0.40 0.80 

TIME 


1.20 1.60 2.00 ms 


Figure 7.19 Simulation of the DET system for the 
linear load changes in Fig. 7. 18 













7.5 SIMULATIONS OF SOLAR-ARRAY SWITCHING SYSTEM 
AND PARTIAL -SHUNT SYSTEM 

A solar array switching system with a buck converter is 
simulated in the shunt mode where the solar array switching 
actually occurs. 

The solar array switching system [NASA2 .MOD] (Fig. 7 . 21 ) is 
configured as shown in Fig. 7. 22. The solar array switching 
unit [AS] is connected with a shunt regulator [FS] and an 
input filter [FI] for the buck converter. The buck converter 
control loop consists of a compensator model [ZP] and a PWM 
model [WM]. Finally, a load model [LO] is connected to the 
buck converter power stage model [BC]. 

The simulation results (Fig. 7.23) begin with a steady 
state in the shunt mode. Since the shunt current is less than 
5A, the solar array switching does not occur during the first 
.2 msec. At 2 msec, the load is step changed so that the load 
current drops sharply. Since the load-power demand is re- 
duced the shunt current increases accordingly. As the shunt 
current reaches the preset value of 5A, the solar array 
switching occurs according to the logic shown in Fig. 4. in 
section 4.2. The responses of bus voltage and buck converter 
output voltage are also given in Fig. 7.23. 

The input requirements data list and the analysis program 
[NASA2.ANC] used for this simulation are given in 
Figs . 7 . 24-25 . 
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Fig. 7.22 Schematic Diagram of [ NASA2 . MOD ] 





0.00 1.00 2.00 3.00 4.00 5.00 6.00 


TIME Ims] 


Fig. 7. 23 Simulation Results of SASU System with Buck Converter 
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INPUT DATA REQUIREMENTS LIST 
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Fig. 7. 24 Input Data Requirements List 


^ NA5A(L AUGl- GIG^NAl ) 

-K- * -fr #■«■* -fi- •« K- i!- -fr ft * i!- -ft *+><-**•»■*** ^^ «• * 

*T1TLE =- SOI A!7 ARHAY SWITCHING SYSTEM WITH ft 

* SHUNT Rl-GULATOH AND A DUCK CONVERTER LOAD ft 

ftftftftftftftftft*ftftftftftftftft#ftftftftftftftft»ftftftftftftftftftftftftftftftft ftftftftftftft 

PARAMETER VALUEG 

R AS « . 001. C AS ^ 5E-5. L AS = IE -6 

Cl AS = 1. E-3> R1 AS 1.0 

NP AS = 32^, NS AS = S8, NNPAS =- 20, DL AS = . 001 

ftftft K■ftft■ftft*ftftftftftftftft^H^•ftftftftft 

INITIAL CONDITIONS 

VB AS = 28. 155. IL AS = 29.479, VCIAG = 28. 184 

ftftftftftftftft 

PARAMETER VALUES 

VR FS = 28. 14, VSAFS = 15. 

R1 FS « 12400, R2 FS = 2. 2E<t>. R3 FS = 1. 784 19 

GM FS 4. VTHFS = 3, ILMFS “ 45 

Cl MC=1,C2 MC=1,C3 MC=1.C4 MC^O 

ft 

INITIAL CONDITIONS 

II FI == 27.985, V2 FI 28. 

PARAMETER VALUES 

R FI=0. 05, L FI=--3E-6. C FI= lOOOE -6 

ftftBC (POWER STAGE! 

INITIAL CONDITIONS 
IL BC=^9. 47, VC DC=20 
PARAMETER VALUES 
L BC~1E-4,C BC=r4E-4 
RL BC=5E--2, RC nC=^ lE -1 

ftftVTG LOOP COM,"EWSATER 
K ZP-. 3, ER 21'= ^. 

WM ZP=4E3 

WZ1ZP='5E4, WZ2ZP-1. 25E2 
WP ZP=1E6 

ft* PWM 

TI WM=20E-6, VP WM^ 6, VQ WM=^. 5 
ER WM=6, CICWM=^0. SCMWM=^0 

ft* LOAD 

RA LO^O. 65, RB LO' . 6. TC L0= lOOE -3 
L L0=lE-6, C L0=^1E~6 
INITIAL CONDITIONS 
XI ZP=34400, X2 2P- . 73.VE ZP=. 3 
IL LO^IO, VC L0= 20 

PRINTER PLOTS 
ONL INE PLOTS 
I NT M0DE=^4 

TMAX=3E-3, TINC==2E 7 
PRATE= 100, 0UTRATE=^100 
SIMULATE 

XTC-X 

PARAMETER VALUES 

RA LO =- 0. 65, RB LO = . 76, TC LO = 2E -3 

DISPLAYl 

V2 BC 

I NT M0DE==4 

DISPLAY5 

VB AS 

DISPLAY2 (OVERPLDT) 

IH FS. 10 AS, II FI 
DISPLAYS, NPPAS 
TMAX=-6E -3, T 1NC=- 2E ~7 
PRATE=^ 50, OUTRATE.^ 100 
SIMULATE 


Fig. 7.25 User' s Parameter Values and Analysis Commands [NASA2.ANC] 
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For the partial shunt system simulation, a similar system 
[PSS6.M0D] (Fig. 7. 26) replaced the SASU model [AS] with the 
partial shunt model [AP] in the previous system configuration 
is developed as shown in Fig. 7. 27. 

The simulation results (Fig. 7.28) begins with a steady 
state in the shunt mode so that the currents and voltages are 
constant for the first 2 msec. Since the load is step changed 
at 2 msec, the load current drops sharply. An increase of 
shunt current results in a decrease of the output voltage 
across the lower portion of the solar array. The voltage 
across the upper solar array increases due to a decrease in 
the load current. As a result, the bus voltage is regulated 
within the specified range. The output voltage of the buck- 
converter reaches the steady- state value after a short 
transient. The input requirements list and the analysis 
program [PSS6.ANC] used for this simulation are given in 
Figs . 7 . 29-30 . 
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PARTIAL SHUNT SYSTEM WITH BUCK CONVERTER LOAD 
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Fig. 7. 27 Schematic Diagram of [PSS6.MOD] 
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with Buck Converter 
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Fig. 7.29 Input Data Requirements List 


ORIGINAL PAGE IS 
OF POOR QUALITY 


TITLE=NASA<LARGE SIGNAL) 


*TITLE = PARTIAL SHUNT SYSTEM WITH BUCK CONVERTER « 
* LOAD « 


PARAMETER VALUES 

R AP = .001, C AP = 500E-6. L AP = lE-6, NP AP = 324 
NSIAP = 29, NS2AP = 29, LL AP =1. , ILMAP = 45 
VIHAP = 15.976, VILAP = 12.57, V2HAP = 15.976, V2LAP = 12.57 
*«■##«•«•«••»■•»•«••«■##«•*•»(■•»!■ *«•*** ■«■ 

INITIAL CONDITIONS 
VD AP = 2S. 161 
ILIAP = 28. 409 


PARAMETER VALUES 

VR FS = 28. 14, VSAFS =15. 

R1 FS = 12400, R2 FS = 2. 2E6. R3 FS = 1.78419 

GM FS = 4. , VTHFS = 3. , ILMFS = 45. 

•fr 


INITIAL CONDITIONS 

II FI = 27.985, V2 FI = 28. 

PARAMETER VALUES 

R FI=0. 05, L FI=3E-6, C FI = 1000E-6 


*JfBC (POWER STAGE) 

INITIAL CONDITIONS 
IL BC=10, VC BC- 20 
PARAMETER VALUES 
L BC=1E"4, C BC=400E--6 
RL 0C=5E-2, RC BC-- lE-1 

COMPENSATOR 
K ZP = O. 3, ER ZP =6 
WM ZP =4E3 

WZIZP =5E4. WZ2ZP = 1 . 25E2 
WP ZP =1E6 


-fi-'K' PWM 

TI WM=20E-6, VP WM=6, VQ WM=. 5 
ER WM=6, CICWM=0. SCMWM=0 

LOAD 

RA LO = 0. 6, RB L0=2. 8, TC L0=100E-3 
L LO =lE-6, C LO=lE-6 
INITIAL CONDITIONS 
IL L0=10, VC L0=20 

I NT M0DE=4 
PRINTER PLOTS 
ONLINE PLOTS 

TMAX=3E-3, TINC=2F-7 
PRATE=10, 0UTRATE=100 
SIMULATE 


XIC-X 

PARAMETER VALUES 

RA LO = 0. 6, RB LO = O. 7, TC LO = 2E-3 
I NT M0DE=4 
DISPLAYl (OVERPLOT) 

VB AP 

DISPLAY2 (OVERPLOT) 

12 AP, IH FS, II FI 
DISPLAYS (OVERPLOT) 

VOIAP, V02AP 
DISPLAY4 
V2 BC 

TMAX=4E -3, TINC=2E -7 
PRATE=100, 0UTRATE=50 
SIMULATE 


Fig. 7. 30 User's Parameter Values and Analysis Commands [PSS6.ANC] 
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7.6 CONCLUSIONS 


In this chapter, various solar array system modes of op- 
eration in the large-signal sense are analyzed using the 
state-plane method. The analysis shows that any equilibrium 
point in Region #2 of Fig. 7.1 is unstable with real 
eigenvalues. The state-plane trajectories are separated by 
a separatrix that passes through Region #2 . The state- 
trajectories converge to either one of the two stable equi- 
librium points located in Region #1 or #3 . Thus the behavior 
of the system is quite different on either sides of the 
separatrix. Since only the equilibrium point in Region #3 
is the desirable operating point, a control method to move 
the operating point from Region #1 to Region #3 is discussed. 
The analysis also suggests a means that leads the system to 
the desired mode of operation under start-up. 

The behavior of the operating point near the solar array 
maximum-power point is analyzed. As far as stability and 
transient response are concerned, the equilibrium points near 
the maximum-power point, in general, are less commendable 
than in the region where the magnitude of the solar array 
output resistance is small, or in the range of a stiffer 
voltage source. 

The solar array system modes of operation including the 
shunt regulator, battery and switching regulator load are 
described. The entire DET system is simulated for a contin- 
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uously varying illumination level and the behavior of the 
system in various modes of operation are analyzed. 
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